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ABSTRACT 



We present the results of a 7 mm spectral survey of molecular absorption lines originating in the disk of a z=0.89 spiral 
galaxy located in front of the quasar PKS 1830—211. Our survey was performed with the Australia Telescope Compact 
Array and covers the frequency interval 30-50 GHz, corresponding to the rest-frame frequency interval 57-94 GHz. A 
total of 28 different species, plus 8 isotopic variants, were detected toward the south-west absorption region, located 
about 2 kpc from the center of the 2=0.89 galaxy, which therefore has the largest number of detected molecular species 
of any extragalactic object so far. The results of our rotation diagram analysis show that the rotation temperatures are 
close to the cosmic microwave background temperature of 5.14 K that we expect to measure at z=0.89, whereas the 
kinetic temperature is one order of magnitude higher, indicating that the gas is subthermally excited. The molecular 
fractional abundances are found to be in-between those in typical Galactic diffuse and translucent clouds, and clearly 
deviate from those observed in the dark cloud TMC 1 or in the Galactic center giant molecular cloud Sgr B2. The 
isotopic ratios of carbon, nitrogen, oxygen, and silicon deviate significantly from the solar values, which can be linked 
to the young age of the z=0.89 galaxy and a release of nucleosynthesis products dominated by massive stars. Toward 
the north-east absorption region, where the extinction and column density of gas is roughly one order of magnitude 
lower than toward the SW absorption region, only a handful of molecules are detected. Their relative abundances are 
comparable to those in Galactic diffuse clouds. We also report the discovery of several new absorption components, 
with velocities spanning between —300 and +170 kms . Finally, the line centroids of several species (e.g., CH3OH, 
NH3) are found to be significantly offset from the average velocity. If caused by a variation in the proton-to-electron 
mass ratio /1 with redshift, these offsets yield an upper limit |^|<4xl0 -6 , which takes into account the kinematical 
noise produced by the velocity dispersion measured from a large number of molecular species. 

Key words, quasars: absorption lines - quasars: individual: PKS 1830—211 - galaxies: ISM - galaxies: abundances - ISM: 
molecules 

1. Introduction As the sensitivity of millimeter instruments has in- 

creased, a large number of molecules have become observ- 
able ingjJ^des_i_nthe_l^ocal Universe (see e.g. IWang et al.l 

A/r 1Kn , , , , 12004 iMartmet al.l [2001 for NGC4945 and NGC253. re- 

More than 150 molecules have been discovered m space, 1 . . : , — r , -n 1 , , , , 1 • , , .„ 

1; r , A j • , • A • . 1 spcctivcly) and molecules arc stdl detected at high redshifts 

as a result ot targeted investigations or spectral surveys s / ' . , , . . , 
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X 



, 1 ^ i ,. 1 • ,11 i (currently up to z=6.42), though only a few species such 

toward Galactic sources such as circumstcllar envelopes v „„ ' „ TT ^^^ 1 „»,>, 1 Utt n , — nr^nnn 
r , j , / TO( ~, , 1(101C r^j — r — n l ormnl as GO. HON, HCO + , and HNG e.g. W alter et al. 2003; 



Patel et aT1l201 lh . mas sive star-forming clouds (e.g. Sgr B2. , . 

Nummelin et aLfemoh. or cold mo lecular clouds fe.g. TMC- T1 ! ose detections are mostly related to ultra-luminous 

1. Ohishi et allll992Tlkaifu et al.ll200l . The study of the § ala f f T ??* Y , unrepresentative of all galax- 

n , ; t ,, n , . les at these epochs. Moreover, line emission from these high- 

rotational and vibrational spectra of these molecules is par- . . . , , 11 , • 

ticularly useful to probe the physical conditions and chem- * g«f**» * dominated by the warmest and densest regions 

. , , , £ . , , rl t» r , • , of the molecular component, 

ical content of the interstellar gas. Most importantly, the ^ 

large number of molecules found in various environments Observations of molecular absorption lines toward a 

demonstrates the richness and variety of interstellar chem- bright radio continuum source allow rare molecular species 

istry. to be detected even in distant galaxies, where they cannot 
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be observed in emission because of distance dilution. This 
technique can thus be used to investigate the chemistry 
and its complexity at high redshifts and trace the chemi- 
cal enrichment history of the Universe. In addition, high- 
z molecular absorbers can serve as cosmological probes of 
the cosmic microwave background (CMB) temperature and 
have attracted much interest in constraining variations in 
the fu ndamental constants of physics (see e.g.. lHenkel et al] 
2009 and references therein). 

Unfortunately, only a handful of high-redshift radio 
molecular abs orbers [J have been discovered to date (see 
the review bv ICombesI 120081). d espite numerous searches 
'e .g. IWiklind fe Combesl Il996al iDrinkwater et~all Il996t 



sarily totally, coupled with the CMB. Rotation tempera- 
tures derived for several molecules are consistent with the 



value T CM B=5.14 K at z=0.89 



iMenten et aIlll999tlHenkel et all 



Combes fc Wiklindl[l997l: 
2009), as predicted for an 



Murphy et all 120031: ICurran et all 12004 l200d~ see also 
Curran et all [20111) . The molecular absorption originat- 
ing in the z=0.89 and z=0.68 intervening galaxies lo- 
cated in front of the quasars PKS 1830—211 (z=2.5) 
and B 0218+357 (z=0.94) respec tively, were discovered by 
IWiklind fc Combesl (1199,1 Il996bfl and have many similar- 
ities: both abso rbers appear to be (faint) nearly face-on 
spiral galaxies (|Winn et al.1 l2002t lYork et all I2005D : both 
galaxies act as gravitational lenses and split the image of 
their background quasar into two main compact compo- 
nents and an Einste i n ring, as seen at r adio wavelengths 
(|Jauncev et al.lll99ll: iPatnaik et aT1ll993l) : they both har- 
bor large column densities of absorbing molecular gas 
(7V(H2) > 10 22 cm -2 ); and in both cases the radio contin- 
uum of the background quasar is bright enough to permit 
sensitive observations with current millimeter instruments. 
A dozen molecular species have indeed been successfully 
detected in these absorbe rs, as well as absorption by their 
rare isotopologues (cf e.g. , Wiklind fe Combeslll995Lfl996bL 
[19981: IMenten et al.lll999t iMuller et al.ll2006D . 

The multiplicity of lensed continuum images of the 
background quasar makes it possible to explore as many 
different lines of sight through the intervening galaxy. In 
the case of B 0218+357, molecular absorption is detected 
only toward one im age of the quasar (|Menten fe Reidlll996t 
IMuller et aill2007lh Toward PKS 1830-211, molecular ab- 
sorption is detected, however, for both of the bright and 
compact images, located one arcsecond apart on the north- 
east and south-west side, r e spectively, of the E instein ring 
(|Wiklind fe Combesl Il998l : IMuller et all [US). The two 
lines of sight intercept the disk of the z=0.89 galaxy (here- 
after referred to as FG0.89, where FG stands for foreground 
galaxy) on either side of its bulge, at galactocentric dis- 
tances of ~2 kpc (to the SW) and ~4 kpc (to the NE). 
An image of the galaxy, obtained with the Hubble Space 
Telescope in the I band is shown in FigfTJ 

The absorption toward the SW image (hereafter 
FG0.89SW) has a large column of absorbing gas (A^(H 2 ) > 
10 22 cm -2 ) a nd is taken as the reference kms" 1 in velocity 
(z=0.88582. IWiklind fc Combesl [l996bh . Remarkably, the 
full width at zero power of the HCO + (2-l) line reaches 
nearly 100 kms -1 for this component, although the galaxy 
is seen almost face-on. We note the relativ ely high kinetic 
temperature of ^80 K (|Henkel et al 20081) and mode rate 
density of a few times 10 3 cm" 3 (|Henkel et all I2009D to- 
ward FG0.89SW, which imply that the excitation of the 
rotational lines is subthermal and mostly, but not neces- 



1 That is, not including the H2/CO absorbers detected 
from optical/UV spect roscopy, such as those discussed in 
iNoterdaeme et al.1 (|2011h . 



adiabatic expansion of the Universe. 

The absorption toward the NE component (hereafter 
FG0.89NE) is located ^147 kms -1 blueward and is nar- 
rower (FWHM^15 kms -1 ), so that both velocity compo- 
nents can easily be resolved kinematically. Nevertheless, the 
column of molecular gas toward the NE image is about two 
orders of magnitude less than in FG0.89SW, and the corre- 
sponding absorption lines are significantly weaker. Little is 
known about the excitation of and the physical conditions 
toward FG0.89NE. However, base d on the analysis o f time 
variations in the NE absorption, IMuller fc Guelinl ((2008) 
argued that it should be caused by a small number (<5) of 
clouds of size ~1 pc and of density a few 100 cm -3 , hence 
resembling Galactic diffuse clouds. 

Molecular absorption has also been detected in the 
galaxies hosting the radio sources B3 1504+377 at z=0.67 
(in C O, HCO+, HCN, and HNC; IWiklind fc Combesl 
Il996bl ). PKS 1413+135 at z=0.25 (in CO, H CO+, 
HCN, HNC, and CN; IWiklind fc Combesl [l99l . and 
PKS 0132-097 a t z=0.76 (although only in OH; 
iKanekar et al.l [20051 ). However, the weakness of the back- 
ground radio continuum sources, narrower linewidths and 
smaller column densities of absorbing gas than the two ab- 
sorbers toward PKS 1830-211 and B 0218+357, ensures 
that it is difficult to perform a deep inventory of their molec- 
ular constituents with current instruments. 

In this paper, we present the first unbiased spectral- 
line survey toward a high-z molecular absorber, located at 
z=0.89 in front of the quasar PKS 1830—211. Observations, 
conducted with the Australian Compact Array and cov- 
ering the whole 7 mm band from 30 GHz to 50 GHz, 
are presented in §f2J Our results are presented in §121 and 
then discussed in § |4j in terms of isotopic ratios (5 14. 1[) . 
the chemical composition relative to other objects (§ 14. 2[) . 
and constraints on the variation in fundamental constants 
of physics (§E5j). Notes for the different molecular species 
detected in the survey are given in appendix. 

2. Observations and data reduction 

Observations were carried out with the Australian 
Telescope Compact Array (ATCA) over three nights, on 
2009 September 1 and 2 and on 2010 March 17 (see 
Table [T]). The array was in the EW352 and H168 configura- 
tion, respectively, and we used the five inner antennas, dis- 
carding the sixth one located 4 km away as we obtained un- 
satisfactory pointing solutions and data quality. The max- 
imum baseline length was then up to 352 m and 192 m 
in either configuration, respectively. With this configura- 
tion, PKS 1830—211 was mostly unresolved, and is further 
considered as a point-like source. The weather was clear in 
both runs. 

For each tuning, the bright quasar 1923—293 was ob- 
served for several minutes to help us perform bandpass 
calibration, and the pointing was checked and updated 
against PKS 1830—211, before integrating. We did not ob- 
serve gain calibrators, as PKS 1830—211 visibilities were 
self-calibrated and the continuum level normalized to unity. 

The Compact Array Broadband Backend (CABB) sys- 
tem was configured to provide a 1 MHz spectral resolution 



2 



Muller et al. 2011: Molecules at z=0.89 




Fig. 1. Image of the z=0.89 galaxy obtai ned with the 
Hubb le Space Telescope in the I band (see also lMuller et al.l 
l2006t ) . The positions of the SW and NE images of the quasar 
are indicated. The distance between them is one arcsec- 
ond. Note the high extinction toward the SW image of the 
quasar, produced by the interstellar medium of the inter- 
vening galaxy. 



per 2 GHz window with dual polarization. Two windows 
were used simultaneously, with central frequencies set as 
reported in Table [Q The whole 7 mm band (30-50 GHz) 
could then be covered in a fairly small number of individ- 
ual tunings. Each part of the band was covered in at least 
two tunings, in order to prevent potential defectuous spec- 
tral channels. Some parts of the spectral band were indeed 
affected by a series of bad channels, as can be seen for ex- 
ample in the final spectrum between 40 and 40.2 GHz in 
Fig lA.ll Some frequency intervals were observed in more 
tunings, to increase the sensitivity, for example at the fre- 
quencies of the deuterated species DCO + and DCN. The 
spectral resolution of 1 MHz corresponds to ^10 kms -1 at 
30 GHz and -6 kms" 1 at 50 GHz. 

The data reduction was performed in two steps. We 
used the MIRIAD package for a first bandpass calibration 
adopting the solution derived from 1921-293 and the self- 
calibration of the PKS 1830-211 data. At this step, the 
analysis of each 2- GHz spectral window revealed a residual 
bandpass, with amplitudes roughly varying by an order of 
one percent over several 100 MHz. To remove this bandpass 
drift residuals - the final sensitivity being of about a few 
times 10~ 3 only -, data were then exported to FITS format 
and subsequently processed with a customised ITT/IDL 
pipeline. 

The pipeline is based on an iterative boxcar drift re- 
moval. Each iteration was designed to construct a model 
spectrum (i.e., a spectrum with any slow variations in the 
bandpass removed) from the raw data following several 
steps: 

— first, for each observed tuning, each polarization and 
each baseline, we flagged outliers in time by assuming 
that the visibilities were constant on timescales longer 
than ten seconds; 

— a drift component was then computed and removed, in 
both time and frequency with different windows widths, 
taking proper care of any flagged or missing data; 

— a first crude spectrum, after evaluating the time- 
averaged visibilities, was obtained and its residuals were 



used to find frequency channels with large deviations 
from a random Gaussian distribution; 

— at this stage, we also examined the spectra for each 
antenna, to check for possible antenna-based problems; 

— the standard deviation of each baseline was then used to 
compute a weighted mean spectrum and the uncertainty 
for each tuning and polarization; 

— finally, a Doppler track correction was applied accord- 
ing to the observatory velocities given in Table [TJ ob- 
tained from the ATCA frequency calculator 0. All ve- 
locities in this paper refer to the Solar System barycen- 
ter rest-frame. We note that the Doppler correction is 
nearly constant during the short observing intervals for 
individual tunings. Data were accordingly re-gridded 
onto a common frequency scale using a one-dimens ional 
(ID) drizzling algorithm (jFruchter fc Hookl[2002h . We 
tested different frequency resolutions and chose a chan- 
nel width of 0.5 MHz. A weighted mean was finally used 
to combine all spectra. 

A model of the spectrum was derived by finding all fre- 
quency channels above a given threshold. This model was 
then removed from the data, and the whole process was 
iterated on the residual datasets. 

We started with a large width for the bandpass resid- 
ual filtering, corresponding roughly to 1 GHz, and a high 
threshold of 10 a to derive the model. The pipeline was run 
until the variance of the resulting spectrum was found to 
have reached a minimum. The model threshold was lowered 
iteratively in steps of two down to 4 a in order to correct 
for the negative corrections of the boxcar drift removal on 
strong lines. Following the same scheme, the filtering width 
was successively divided by two until it reached a width of 
16 MHz. 

This whole procedure ensures that the second order 
bandpass is removed, without making any a priori assump- 
tions about the presence of strong absorption lines. The 
resulting spectrum is shown in Fig lA.ll 

3. Results and analysis 

An overview of the 20-GHz-wide resulting spectrum is 
shown in Fig [2] The frequency coverage corresponds, in the 
rest frame of the 2=0.89 galaxy, to the range 57-94 GHz. 
The lower end of this frequency interval is relatively unex- 
plored because of poor atmospheric transmission caused by 
a forest of fine-structure lines of O2 around 60 GHz. The 
rms noise level is not uniform over the whole band, and 
varies owing to the higher system temperatures at higher 
frequencies and the different integration times. The full sur- 
vey is presented in Fig lA.ll 

Only a handful of absorption lines reach an optical 
depth of 10% (see EqfT|) or more. They originate in the 
SW absorption component and correspond, in order of 
decreasing optical depth, to: HCO+, HCN, HNC, C 2 H, 
N 2 H+, H 2 CO, c-C 3 H 2 , and H 13 CO+. A ll these species have 
been detected in previous observations ()Wiklind fe Combes! 
119981: iMenten et all 1 19991 ) . The weaker NE absorption at 

V 147 km" 1 is detected only for HCO+, HCN, HNC, 

C 2 H, c-C 3 H 2 , and H 2 CO. 

The major result of this survey, though, is the detection 
of a collection of weak absorption lines toward the SW com- 

2 http:/ /www. narrabri.atnf.csiro.au/observing/obstools/velo. html 
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Table 1. Journal of the observations. 



Freq. 


Time 


At 


Vobs 


(GHz) 


(UT) 


00 


(km s" 1 ) 


2009 September 1 


31 & 33 


6:49 


0.8 


+25.71 


37 & 39 


8:14 


1.2 


+25.84 


44 & 47.5 


10:00 


1.5 


+26.03 


42 & 46.5 


11:56 


1.6 


+26.23 


35 & 38 


13:32 


0.9 


+26.36 


2009 September 2 


32 & 34 


6:44 


0.8 


+25.95 


41 & 45 


8:08 


1.2 


+26.07 


46.5 & 49 


10:19 


1.4 


+26.30 


43 & 47 


11:55 


1.0 


+26.47 


36 & 38.5 


13:10 


0.7 


+26.57 


2010 March 17 


37.25 & 40 


19:16 


3.0 


-29.44 


46.1 & 48.1 


22:22 


3.0 


-29.16 



Notes. Frequencies indicate the center of each 2 GHz 
ATCA/CABB bands; time is given at the middle of the time 
interval At spent on the frequency tuning, and V b s is the obser- 
vatory velocity used to perform Doppler correction (with respect 
to the Solar System barycenter rest-frame). 



were also detected. Elemental isotopic ratios of C, N, O, 
and Si are derived and discussed in £14.11 
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Fig. 3. Cumulative number of lines detected toward 
FG0.89SW with optical depths higher than a specified 
value. 



ponent, down to a few times 10 -3 of the continuum inten- 
sity, as illustrated in FigGH where the cumulative number 
of detected lines is plotted as a function of opacity. These 
lines were identified using the on-lin e Cologne Databas e 
for Molecular Spectroscopy (CDMS. iMuller et al.1 [20(51 . 
and Jet Propulsi on Laboratory Mo lecular Spectroscopy 
Database (JPL. iPickett et al.ifl998t ). The corresponding 
transitions are listed in Table IA.1I 

All detected transitions (except the K > 1 transitions 
of CH3CN, see § IB. 141) are ground state or have low energy 
levels. Most molecules have high dipole moments, and the 
gas density is too moderate for collisions to play an efficient 
role in the excitation. Rotational transitions are therefore 
nearly in radiative equilibrium with the CMB, which has 
a temperature of ^5 K at z=0.89. Several species are de- 
tected in more than one transition, and a rotation diagram 
analysis indeed yields rotation temperatures in agreement 
with this value. A low energy level is therefore a good filter 
to identify transitions/species in the molecular databases, 
and this is why we believe that our identifications are ro- 
bust, even for species with a single detected transitions. The 
confusion is also low, with a rough average of four clear lines 
per GHz over the 20-GHz-wide band. 

The census of molecules detected to date in the z=0.89 
galaxy located in front of PKS 1830—211 is given in Tabled 
It is quite remarkable that, despite its distance, this z=0.89 
galaxy has now become the object with the largest num- 
ber of detected molecules outside the Milky Way. A large 
fraction of the total number of these species (28 out of 34) 
are detected in our survey, among them 19 for the first 
time toward this source. As far as we know, this work also 
represents the first extragalactic detection for SO + , I-C3H, 
1-C 3 H 2 , H 2 CCN, H 2 CCO, C 4 H, CH 3 NH 2 , and CH 3 CHO. 
In addition, isotopologues of HCO+, HCN, HNC, and SiO 



3 http:/ /www.astro. uni-koeln.de/cdms 

4 http:/ /spec.jpl. nasa.gov/home. html 



3.1. Fitting procedure 

We fitted a Gaussian profile to all the transitions of a given 
species with a single velocity and a single linewidth. The 
ortho and para forms, when relevant, were considered as 
the same species for the fit. The relative intensities of all 
species with resolved or marginally resolved hyperfine struc- 
ture were fixed, as it was observed that, in the case of 
C 2 H and HCO (for which the hyperfine structure is clearly 
resolved), they follow the predicted ratios assuming local 
thermal equilibrium (LTE). The velocity and linewidth of 
isotopic variants were fixed relative to the main isotopo- 
logue, except for HCO + and HCN, the lines of which are 
likely saturated. For these, the corresponding isotopic vari- 
ants were instead fixed relative to H 13 CO + and H 13 CN. 
For some species with data for the various transitions of 
low signal-to-noise ratio (S/N), we decided to fix the ve- 
locity (e.g. 1-C3H 2 ) or the linewidth (e.g. SO + ) after initial 
free fits. Finally, the fit for CH3CN is degenerate because of 
the blending of the K=0 and K—l transitions, and for this 
species, we therefore chose to adjust a synthetic spectrum 
by eye. The fitting results are given in Table IA.1I 

3.2. Time variations 

A monitoring of the HCO + J=2-l absorption profile 
betwe en 1995 and 2008 rev ealed significant time varia- 
tions (jMuller fc Gu clin 2008|) , in particular with the quasi- 
disappearance of the —147 kms -1 absorption component 
between 2006 to 2007. In the HCO+ J=2-l spectrum 
acquired by Muller & Guclin in October 2008 with the 
Plateau de Bure interferometer, this component had re- 
turned with an absorption of ~15% of the total continuum 
intensity. Variations were found to be correlated between 
the —147 kms -1 component and the blue wings of the 
kms -1 absorption, hence should originate from morpholog- 
ical changes in the background quasar. 

From Very Long Baseline Interferometry continuum 
measurements, the relative distance between the NE and 
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SW lensed images of the quasar was found to vary by 
as much as 20 micro-arcseconds within a few months 
(|Jin et al.l 120031 ) . iNair et all (|2005l) proposed that this ap- 
parent motion is caused by the sporadic ejection of plas- 
mons along a helical jet by the quasar. Consequently, 
the continuum illumination scans different pencil beams 
through the foreground absorbing clouds, where the appar- 
ent displacement of the quasar lensed images corresponds 
to a projected distance of about a parcsec, i.e. a scale com- 
parable to the size of Galactic molecular-cloud cores. 

The relative differences (normalized to the averaged 
spectrum) for the spectral regions covered in both runs 
(September 2009 and March 2010) are shown together with 
the averaged spectrum in Fig lA.ll They are at most ~5% 
for the HCO + , HCN, and HNC lines, and barely visible for 
all other lines. Given the small relative changes, data from 
both observing sessions were combined to improve the S/N 
of the data. 



3.3. Background continuum illumination 

Absorption intensities (1^), measured from the continuum 
level, were converted to optical depths (t„) according to 



= - In 1 - 



fJi 



hi 



(1) 



where I^g is the intensity of the background contin- 
uum source, and f c the source covering factor. At ra- 
dio wavelengths, the continuum emission of the quasar 
PKS 1830—211, lensed by the z=0.89 intervening galaxy, 
is dominated by two bright and compact components (to 
the NE and SW), separated by 1" and included in an 
overall fainter Einstein ring. The flux ratio of the NE to 
SW compact c omponents is ~1.5-1.6 at cm wavelengths 
as measured bv ISubrahmanvan et al.l (jl990) . The remain- 
ing emission (Einstein ring and fainter other components) 
contributes only to a few percent of the total continuum 
intensity. 

Molecular absorption in the z=0.89 galaxy is observed 
toward both the SW and NE continuum images. The HCO + 
J=2-l absorption toward the SW image shows a fiat- 
bottom profile, indicating saturation. Assuming that, at 
this wavelength, all the continuum emission comes from 
the NE and SW images of the quasar, it is then possible 
to derive the relative flux NE/SW of both images from the 
level of this saturation, without actually resolving them. 
The NE/SW ratio was measured to be about 1.7 with a 
dispersion of 0.3 over several ob servations between 1995 
and 2007 dMuller fc Guelinl l2008h . i.e. roughly similar to 
the flux ratio directly measured at cm wavelengths. Some 
of the fluctuations might be due to morphological changes 
in the quasar (plasmon burst, microlensing) and the time 
delay between the two line of sights. 

Our 7 mm survey is limited in terms of spectral resolu- 
tion, and the flat-bottom part of the HCO+/HCN J=l-0 
lines (supposing that they are also saturated) is smooth, 
preventing the same analysis. Nevertheless, as the HCO + 
1 — absorption reaches an apparent depth of ^35% 
and 8% toward the SW and NE components, relative to 
the total continuum level, we can certainly state that 
0.35 < Ib s (SW) < 0.92 and 0.08 < I hg (NE) < 0.65. 
If, in addition, we impose the ratio R =NE/SW, we get 
the stronger constraints of 0.35 < I bg (SW) < ^ 



0.35xi? < 7 6s (NE) < 0.65, which gives I b9 (NE)=0.61±0.05 
and 7 fcff (SW)=0.38±0.03 with #=1.6. We adopt these val- 
ues. The rough uncertainty in the illuminating background 
continuum for both NE and SW components is therefore 
less than 10%. The residual contribution from components 
other than the NE and SW images is de facto less than 
[1 - 0.35 x (1 + R)], that is 9% with R=1.6. 

This is valid, of course, if no changes happened in the 
quasar within the time corresponding to the time delay be- 
tween the two lines of sight (^24 days) before our observa- 
tions. Unfortunately, we do not have flux monitoring data 
corresponding to these epochs. Nevertheless, that we do not 
see significant variations in the absorption line profiles be- 
tween the observations of September 2009 and March 2010 
is reassuring. 

We can now investigate how this uncertainty in the 
background continuum illumination propagates to the 
derivation of opacities, by means of Eq|T] In the case of 
weak lines (I « Ib g ), we get = -7-^, and our opaci- 
ties are therefore estimated more accurately than 10%. For 
the other lines, we get = ^j 3 - x F(I/Ib g ), where the 
function F(x) = x/((x — l)ln(l — x)), and tends to infin- 
ity for x — > 1, although the rise is slow, as F(0.9) (that is 
I/hg = 90%) is still < 4. For example, the corresponding 
uncertainty in HNC opacity (the third strongest absorption 
line in our survey, after HCO + and HCN) is less than 20%. 

Finally, we assume here (as in the rest of the paper) a 
source covering factor f c of unity for all absorption compo- 
nents. It is however possible that f c < 1 and changes with 
frequency, as the continuum emission probably becomes 
more e xtended at lower f requencies (see also the discus- 
sion bv lHenkel et al.| [2009). As a result, the optical depths 
derived from Eq[T] are, strictly speaking, lower limits. We 
note however that the HCO + and HCN J=2-l lines toward 
the SW absorption are optically thick and saturated. The 
source covering factor for this component is thus certainly 
close to unity. Toward the NE absorption, the situation is 
not so clear. The formal observational constraints from the 
HCO+ 1-0 line are r > 0.14 and f c > 0.13. While the non- 
detection of H 13 CO+ (with HCO+/H 13 CO+>71) suggests 
that the line is likely optically thin, the source covering 
factor is still only poorly constrained. 



3.4. LTE analysis and rotation diagrams 

For optically thin lines, and assuming a Boltzmann distri- 
bution characterized by a rotation temperature T rot , col- 
umn densities N co i can be derived as 



3h Q{T ro t) exp ( 



k B T T , 



LTE 



JttVSj [i_ exp (^-)] 



rdV, 



(2) 



where Q(T rot ) is the partition function, Ei the energy of 
the lower level with respect to ground state, /i the dipole 
moment, v the frequency, and S u i the line strength. Under 
these conditions, the column density of a detected molecule 
can be calculated directly with the knowledge of only one 
physical parameter: the rotation temperature T rot . 

Since some species are detected in several transitions, an 
excitation analysis can be done. Under the Rayleigh- Jeans 
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(RJ) approximation {hv << fc^T), EqJ5]canbe re-arranged 
as 

log ( 3k B jTdV \ = ( N col \ _ mogje) 
\8ir 3 vn 2 S u i) \Q{T rot )T rot J k B T rot 

Plotting the quantity on the left side of this equation ver- 
sus the energy of the lower level for each observed transi- 
tion then allows the rotation temperature T rot to be mea- 
sured as the inverse of the slope of a linear fit to the data 
points in the rotation diagram. Column densities were also 
derived, following Eq|21 from the intercept of the y-axis 
in the rotation diagram. Because the uncertainties in the 
data points in the rotation diagrams do not follow a normal 
distribution, the rotation temperatures, column densities, 
and their corresponding uncertainties were estimated using 
Monte Carlo simulations. We note that the RJ approxima- 
tion is still valid for T > 4.5 K at 94 GHz, the highest 
rest-frame frequency reached in our survey. 

The rotation diagrams for molecules detected toward 
FG0.89SW in several transitions are presented in Fig|4j and 
the resulting measurements of rotation temperatures and 
column densities are listed in Table 2] These rotation tem- 
peratures are plotted in Figj5] as a function of the molecu- 
lar dipole moment. The results for column densities toward 
FG0.89NE are given in Table 

Most of these species have high dipole moments, and 
consequently high critical densities (>10 5 cm -3 ). For low 
density gas (~10 3 cm -3 ), the population of the different 
energy levels is not efficiently thermalized by collisions, but 
is coupled with the ambient radiation field, which, in the 
absence of any other radiative excitation (UV or IR pump- 
ing), consists of CMB photons. In this case, rotation tem- 
peratures then nearly equal the CMB temperature, which is 
expected to be Tcmb=5.14 K at z=0.89, assuming a scal- 
ing proportional to (1 + z). The derived rotation tempera- 
tures are consistent with this value. A thorough excitation 
analysis is beyond the scope of this paper, as it would re- 
quire a careful treatment of collisional rates and partners 
(H2 and e - ) for each molecule. This will be the subject of 
a future paper, in which we will also study the excitation 
in FG0.89NE with new multi-transitions observations, in 
order to derive an accurate and robust measurement of the 
CMB temperature at z=0.89 through the two independent 
lines of sight. 

As a first step to interpreting the observations, we pro- 
duced a LTE synthetic spectrum resulting from the combi- 
nation of absorption lines from the molecular species listed 
in Table 2] We assume a common rotation temperature 
Irot=Tc7MB=5.14 K for all species and for the sake of sim- 
plicity, the absorption profile is set to be Gaussian with a 
FWHM of 20 kms" 1 centered at V=0 kms -1 for the SW 
absorption (some lines are slightly offset from these values, 
see Table Hand 34~3)) and at V=-147 kms -1 for the few 
species for which NE absorption is detected. The partition 
functions were calculated for T=5.14 K by performing a 
direct sum over all energy levels. Their values are given in 
Table Q] Following EqJ^l the only free parameter is then 
the column density of each species. A kinetic temperature 
of 50 K was adopted to take into account the ^-ladder pop- 
ulation of the symmetric-top molecule CH3CN (see i|B.14[) . 

The resulting LTE spectrum, obtained with the LTE 
column densities listed in Table [H is overlaid on top of the 
observed spectrum in Fig lA.ll Despite the simplicity of the 
model, it reproduces relatively well the observations. 
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Fig. 5. Rotation temperatures as a function of molecu- 
lar dipole moment. The dashed line indicates the value 
T'cmb=5.14 K at z=0.89, predicted from standard cosmol- 
ogy. Error bars indicate the 95% confidence interval. 



In addition, we list in Table [5] the upper limits to the 
column densities and abundances obtained for some in- 
terstellar species of interest that have low-energy transi- 
tions within the frequency coverage of our survey. The up- 
per limits to the integrated intensities were calculated as 
3<t t ^/<5V AV1/2, where a r is the optical depth (EqllJ noise 
level, 5V is the velocity resolution, and AV^ the linewidth 
at half-maximum, which we fixed to 20 kms -1 . The upper 
limits to the column densities were calculated from Eqj2l 
again assuming LTE and T rot =5.14 K. 

3.5. Additional velocity components 

Previous observations have shown that the molecular ab- 
sorption profile consists of two main velocity components, 
separated by 147 kms -1 and located in front of the SW 
and NE images of the quasar (see Introduction). However, 
this picture should now be updated with the discovery of 
additional absorption components. Close to the redshifted 
frequencies of both the HCN and HCO + 1-0 lines (red- 
shifted to ~47 and 47.3 GHz respectively) in Fig lA.ll sev- 
eral additional weak lines can be seen, which we failed to 
identify with other species. Those lines clearly arise at the 
same velocities for both HCN and HCO + 1-0 transitions, 
as shown in Fig|6l They therefore correspond to additional 
velocity components of HCN and HCO + 1-0 in the z=0.89 
galaxy. The new components are located at —300, —224, 
—60, and +170 kms -1 (in addition to the previously known 
ones at and —147 kms -1 ). They probably arise in front 
of secondary continuu m peaks along the Ein stein ring (see 
for example Fig. 3 bv IChengalur et all [1999) . It would be 
interesting to determine their locations to constrain the 
kinematics of the galaxy, as well as their background con- 
tinuum illumination to derive the associated opacities and 
column densities of gas along different sightlines through 
the disk. This is unfortunately impossible with the current 
data, which are limited in terms of angular resolution and 
sensitivity. Interestingly, the —300 and +170 kms -1 com- 
ponents ha ve velocities beyond the range covered by HI 
absor ption (jChengalur et al.lll999t iKoopmans fc de Bruvnl 
|2005[ ) , and are nearly symmetrical with respect to the cen- 
tral velocity (^—80 kms -1 from the velocity of the SW 
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absorption). These two components may arise from rapidly 
rotating gas near the center of the galaxy, e.g., in a circum- 
nuclear ring. 

Finally, we emphasize that we were unable to identify 
any lines o riginating from th e second intervening galaxy at 
z=0.1926 (jLovell et al.lll996|). neither from the host galaxy 
of the quasar at z=2.5 ( Lidman et al.l Il999t ). nor from 
Galactic absorption. 
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Fig. 6. Spectra of the HCO+ and HCN 1-0 lines, showing 
the absorption components at and —147 kms -1 (previ- 
ously known, and arising toward the SW and NE images 
of the quasar, respectively), as well as the additional ab- 
sorption components at velocities —300, —224, —60, and 
+170 kms -1 . 



4. Discussion 

The advantages of our survey are manifold. First, the 4-mm 
rest- frame window (redshifted to 7 mm at z=0.89), includes 
numerous ground state and low energy level transitions 
of a large range of intermediate-size interstellar molecules, 
and is therefore particularly interesting for an absorption 
study. Lighter molecules, e.g. hydrides, on the one hand, 
and heavier and more complex molecules (with > 6 atoms) 
on the other hand, generally have these transitions at higher 
(mm/submm window) and lower (~cm wavelengths) fre- 
quencies, respectively. The achieved sensitivity, of a frac- 
tion of a percent relative to the continuum intensity, and the 
large spectral coverage (20 GHz band) then allows us to de- 
tect many weak lines of rare molecules (i.e. with abundances 
of a few 10 -11 relative to H2) and in several cases, to build 
their rotation diagrams. For both runs (September 2009 
and March 2010), observations were carried out within one 
or two days, i.e. a ti mescale shorter than the time delay of 
the le ns (~24 days. lLovell et alj|1998t IWiklind k Combes! 
Il999h and much smaller than the timescale for variations 
in the absorption profi les, of the order of several months 
(|Muller k Guelip| l2008). The comparison of lines in the fre- 
quency overlap between the two datasets shows that the rel- 
ative variations between September 2009 and March 2010 
are smaller than -5% for the strong HCO+/HCN 1-0 lines 
and yet smaller for other lines. Time variations are therefore 
negligible, and moreover, would affect only the few lines in 
the small frequency overlap between both observing runs. 
Finally, the absorption technique provides a high spatial 



resolution, equivalent to the angular size of the quasar con- 
tinuum images and corresponding to a few parsecs in the 
plane of FG0.89. 

On the other hand, several limitations have to be kept 
in mind. First, the detection of molecular species in our 
survey (for a given abundance) is obviously biased towards 
those that have i) a relatively high dipole moment, ii) a 
relatively low partition function, and Hi) low energy tran- 
sitions within the frequency coverage (57-94 GHz, rest- 
frame). The molecular inventory is then limited by a com- 
bination of these properties. Second, the angular resolution 
of our observations is insufficient to resolve the background 
continuum emission. While the different absorption com- 
ponents are well-resolved kinematically, this introduces an 
uncertainty in the continuum illumination, hence also in 
the line opacities and column densities. We argue however 
that the associated uncertainties are only ~10% for the 
bulk of the lines (see ? !3.3p . Moreover, the spectral resolu- 
tion is only a factor of a few smaller than the width of the 
line (~20 kms -1 ). This reduces the absorption line depths, 
likely causing opacities to be slightly underestimated. We 
expect that this affects weak lines the most. The limited 
spectral resolution also prevents us from measuring the sat- 
uration level of the HCO+ /HCN 1-0 SW lines, which could 
help us to derive the NE/SW flux ratio. Finally, we assume 
a source covering factor of unity for each absorption compo- 
nent, implying that the derived optical depths and column 
densities are, strictly speaking, lower limits. 

4.1. Isotopic ratios 

Isotopic abundances are directly affected by nucleosynthe- 
sis processes in stellar interiors, and isotopic ratios are 
thus good probes of the chemical enrichment history of the 
Universe. For this purpose, molecular absorption offers in- 
teresting prospects, owing to the rotational transitions of 
different isotopologues being easily resolved at mm wave- 
lengths, such that low column densities of the rare isotopo- 
logues are still detectable e.g. for most isotopes of C, N, 
O, and S elements, and that abundances can be inferred 
directly from absorption optical depths (if the line of the 
most abundant isotopologue is not saturated). 

The 12 C/ 13 C, 14 N/ 15 N, 16 0/ 18 0, 18 Q/ 17 Q and 32 S/ 34 S 
isotop ic ratios were derived in FG0.89SW by I Muller et al.l 
(|2006l ). based on observations of the different isotopologues 
of HCO+, HCN, HNC, CS, and H 2 S at 3 mm with the 
Plateau de Bure interferometer. These ratios, particularly 
18 0/ 17 and 32 S/ 34 S, reveal significant differences when 
compared to source s in the Local Universe (see Table 7 from 
iMuller et al.l l2006f) . Interestingly, comparable values are 
found using the same absorption technique in the z=0.68 
galaxy located in front of the quasar B 0218+357 (Muller 
et al. in prep.). These results are essentially consistent with 
the expectation that low-mass stars (<1.5 M ) had no (or 
only a short) time to play a major role in the gas enrichment 
of such young (a few Gyr old) galaxies. 

Our 7 mm survey now allows us to check the previous re- 
sults obtained at 3 mm in FG0.89SW, in particular for lines 
with different optical depths. In addition, time variability 
is not a concern for this 7 mm dataset. Unfortunately, the 
limited spectral resolution of the survey prevents us from 
measuring ratios across the lines, which is a thorough test 
of possible saturation, and we simply derive the isotopic 
ratios from the ratios of integrated opacities. 
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The spectra of different isotopologues of HCO + , HCN, 
HNC, and SiO as observed in our 7 mm survey are shown 
in FigjT] Results for the corresponding isotopologue abun- 
dance ratios are given in Table where we also report the 
ratios previously measured at 3 mm. Our final estimates 
for the various isotopic ratios in FG0.89SW are given in 
Table [3] The values obtained from 7 mm and 3 mm data 
are in very good agreement. We hereafter discuss how we 
derived those ratios. 
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Fig. 7. Spectra of the different isotopologues of HCO + , 
HCN, HNC, and SiO. 



4.1.1. Deuterium 

We do not detect lines from any of the deuterated species 
DCO+/DCN, down to levels of a few 1CT 4 , which is 
higher than the D/ H cosmic abundance of 2.5 xlO -5 (e.g. 
ISpergel et alJl2003» . This is consistent with no strong frac- 
tionation enhancement, as expected in a relatively warm 
gas component (Tfcj n ~5Q K) with a moderate density, n(H 2 ) 
of the order of a few 10 3 cm~ 2 . 

4.1.2. Carbon 

It is not straightforward to go from the abundance ratio 
12 C*/ 13 C*, where C* is a given carbon-bearing molecule, to 
a value of the 12 C/ 13 C isotopic ratio, because the measured 
abundance ratio can be affected by opacity effects, chemi- 



cal fractionation and/or selective photodissociation. Hence, 
observations of several 12 C- and 13 C- bearing molecules are 
highly desirable. 

The interstellar gas phase abundance of 13 CO and 13 C + 
relative to CO and C + , can be affected by the isotopic frac- 
tionation reaction 

13 C + + CO < — > 13 CO + C + + 35 K. (4) 

In cold molecular clouds, 13 C then becomes mostly trapped 
in 13 CO, and the abundance of the 13 C + ion decreases. 
The 12 C*/ 13 C* ratio for species formed from CO should 
thus reflect the 12 C/ 13 C isotopic ratio, while for others 
12 C*/ 13 C* > 12 C/ 13 C should be ob served. This was first 
investigated by lLanger et al.l (|1984[ ) for various carbon- 
bearing species. They found that the 12 C/ 13 C isotopic ra- 
tio should be within the limits set by the 12 CO/ 13 CO and 
H2 2 CO/H3, 3 CO ratios. Since HCO + can be produced from 
both CO and other species, it can be enhanced in either 
12 C or 13 C, depending on the physical conditions. 

The 12 C*/ 13 C* abundance ratios are listed in Table [7J 
for relevant species. The HCO + /H 13 CO + ratio might be 
underestimated because of the large opacity of the 12 C 
isotopologue. The opacity of HNC, on the other hand, is 
lower, and the HNC/H 13 NC abundance ratio might more 
accurately reflect the 12 C/ 13 C elemental ratio. We do not 
detect transitions from either the 13 CCH or C 13 CH iso- 
topomers, which implies that the lower limits to the abun- 
dance ratios are C 2 H/ 13 CCH>58 and C 2 H/C 13 CH>43. 
iSakai et al.l ([2010) found 13 C abundance anomalies in the 
13 C-isotopomer of C 2 H toward TMC 1 and L1527, with 
the 13 C-species underabundant relative to the interstel- 
lar 12 C/ 13 C ratio. In addition, the two carbon atoms do 
not appear to be equivalent in the formation pathways 
of the molecule, as ISakai et al.l (|2010f ) measured an abun- 
dance ratio C 13 CH/ 13 CCH = 1.6 toward both sources. 
They conclude that carbon-chain molecules are not indi- 
cated to determine the 12 C/ 13 C elemental ratio, because of 
the positional differences and heavy 13 C dilution. The non- 
detection of H^ 3 CO points to a ratio H 2 CO/H^ 3 CO>173, 
suggesting that fractionation could affect our estimate of 
the 12 C/ 13 C isotopic ratio. On the basis of all these results, 
we choose to adopt the average ratio 12 C/ 13 C=35 ± 11, 
obtained from the HCO+, HCN, and HNC species. 

This value of 12 C/ 13 C is les s than half the ratio mea- 
sured in the Solar System (89. lLoddersI 120031) . but closer 
to that derived in the local interstellar medium (59 ± 2, 
iLucas feTL iszt 199^) and especially close to th at in the star- 
burst galaxies NGC253 and NGC4945 (40-5 0. iHenkel et all 
Il993t ICurran et al.ll200lHWang et al.ll2004h . It is, however, 
much lower than the value that would be expected in a 
poorly proces sed environment (e.g . iKobavashi et al-lfeoilt ). 
We note that iMartm et all (|201Cft used C 2 H and CO iso- 
topologues to revisit the 12 C/ 13 C ratios in NGC253 and 
M82 and found values larger than the ratios measured in 
previous studies. The interstellar 12 C/ 13 C isotopic ratio 
clearly remains difficult to determine. 

Toward the NE component, we do detect neither the 
H 13 CO + nor the H 13 CN 1-0 absorptions, with lower limits 
HCO+/H 13 CO+>71 and HCN/H 13 CN>41 (at 3 a level), 
respectively. As the NE absorption is located at an even 
larger galactocentric distance than the SW, thus in a region 
likely to be even less processed, we expect a 12 C/ 13 C ratio 
larger than that toward FG0.89SW. 
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4.1.3. Nitrogen 

There are two ways of estimating the 14 N/ 15 N isotopic ra- 
tios from our data: either from the direct HCN/HC 15 N 
and HNC/H 15 NC abundance ratios, or, indirectly from the 
H 13 CN/HC 15 N and HN 13 C/H 15 NC ratios, after correcting 
for the 12 C/ 13 C isotopic ratio. 

We obtain the ratios HCN/HC 15 N=192 ± 11, and 
HNC/H 15 NC=340 ±^°. The former ratio might be un- 
derestimated because of opacity effects on the HCN line, 
but the second ratio has a large uncertainty. Using the 
13 C variants, we obtain H 13 CN/HC 15 N=5.3 ± 0.3 and 
HN 13 C/H ls NC=7.3 2 _f 5 . The use of a large 12 C/ 13 C ra- 
tio (>50) gives a 14 N/ 15 N ratio close to the value observed 
in th e local interstellar medium ( 
1998) and in the Solar System 



=237 
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( 



Using our estimate 12 C/ 3 C=35 



+27 

=272. 2] |Loddersl 12005 
± 11 and the ratio 
H 13 CN/HC 15 N=5.3 ± 0.3, we derive 14 N/ 15 N=190 ± 60. 

Whether this value actually reflects the true 14 N/ 15 N 
isotopic ratio is however uncertain i f fractionation is an 
issue for nitrogen-bearing species (iRodgers fc Charnlevl 
l2008allbh . The recent measurements of iLis et all (|2010D to- 
ward the two cold dense molecular clouds Barnard 1 and 
NGC1333 do not suggest a high 15 N enhancement. 

We note that the 15 N isotopologues of N2H + are not 
detected, implying formally that 14 N/ 15 N > 90. 



4.1.4. Oxygen 

As for 14 N/ 15 N, the 16 0/ 18 isotopic ratio can be esti- 
mated either directly from the HCO + /HC 18 + abundance 
ratio, or indirectly from H 13 CO+/HC 18 0+. The former 
method yields a ratio 16 O/ 18 O=50± 1, while the latter im- 
plies that 16 O/ 18 =80 ± 25. Even considering that our value 
of 12 C/ 13 C could be underestimated, the 16 0/ 18 isotopic 
ratio derived in FG0.89SW is much lower th an that in both 
the l ocal interstellar medium f=6 72 ± 110. jLucas fc Lisztl 
1998) and the Solar System (=499,|Iodders|l2Q03D, in likely 
connection with the youth of the galaxy. 

Remarkably, the 18 0/ 17 isotopic ratio, measured di- 
rectly from the H 18 CO+ and HC 17 0+ J=l-0 optically thin 
lines, also differs significantly from valu es for sources in th e 
Local Universe (< 6, see e.g. Table 7 of lMuller et"aDl2006h . 
While a value of 18 0/ 17 0=12 l" 2 was measured from the 
J=2-l lines, the marginal detection of HC 17 + at 7 mm 
(FigO yields H 18 CO+/HC 17 O+=20 t\- This new estimate 
of the 18 0/ 17 isotopic ratio confirms the large value pre- 
viously obtained from the J=2-l lines at 3 mm. 



(jAgafonova et al.ll201lh . and also suggest that there is a 
significant relative overabundance of heavy Mg isotopes. 
Silicon and sulfur (and magnesium) are both produced by 
massive stars. That the neutron-enriched isotopes appear 
to be more abundant at z=0.89 than in the Local Universe 
might provide some constraint on their nucleosynthesis. 



Table 7. 

FG0.89SW. 



Isotopologues abundance ratios toward 



Abundance ratios 


@7 mm 


@3 


mm 


DCO+/HCO+ 


< 7.1 x 10 -4 






DCN/HCN 


< 8.4 x 10~ 4 






DNC/HNC 




< 


.015 


DCO+/H 13 CO+ 


< 0.015 






DCN/H 13 CN 


< 0.030 






HCO + /H 13 CO + 


21.2 ±0.2 


28 


±3 


HCN/H 13 CN 


36.3 ±0.5 


40 


+7 

-5 


HNC/HN 13 C 


47 ± 2 


27 


±3 


H 2 CO/Hf CO 


> 173 (3a) 






C 2 H/ 13 CCH 


> 58 (3a) 






C 2 H/C 13 CH 


> 43 (3a) 






HCN/HC ib N 
HNC/H 15 NC 


192 ± 11 






340 tlf 


166 


+194 
-58 


H 13 CN/HC 15 N 
HN 13 C/H 15 NC 


5.3 ±0.3 


4.2 ±0.3 


7.3 ±™ 


6 


+4 
-2 


HCO + /HC lts O + 


48 ±1 


53 


+16 

-10 


H 13 CO+/HC 18 0+ 


2.29 ±0.05 


2.01 ±0.07 


HC 18 + /HC 17 0+ 


19.6 ±1 


12 


+3 
-2 


28 SiO/" H SiO 


11 ±t 







Table 8. Isotopic ratios toward FG0.89SW. 



Isotopic ratios 


7 mm 


3 mm 


Averaged 


D/H 


< 7 x 10" 4 


< 0.015 




12 C/ 13 C 


35 ±11 


32 ±6 


33 ±5 


14 N/ 15 N 


190 ± 60 


143 ± 30 


152 ± 27 


16 o/ 18 o 


80 ±25 


66 ± 12 


69 ± 11 


18 0/ 17q 


20 tl 


12 tl 


13 tl 


32gy34g 




10.5 ±0.6 


10.5 ±0.6 


28 Si/ 29 Si 


111 4 




lit 4 



Note s. The values at 3 mm were rederived from iMuller et al.1 
(2006) data using the same methodology as described in the 
text. 



4.1.5. Silicon 

The detection of the 28 SiO and 29 SiO J=2-l line allows us 
to estimate the isotopic ratio 28 Si/ 29 Si=ll I 4 - Both SiO 
lines are optically thin, unless the source covering factor is 
much lower than unity. This 28 Si/ 29 Si ratio at z=0.89 is 
nearly half the value measured in the Solar System, 19.6 
(|Loddersll200"l . 

The situation is similar to that for the 32 S/ 34 S 
isotopic ratio, estimated to be 10.5 ± 0.6 from the 
CS/C 34 S and H 2 S/H^S abundance ratios in FG0.89SW 
(jMuller et al.ll2006h . whereas it is 22 in the Solar System. 
Interestingly, the magnesium isotopic ratios were derived 
from UV-spectroscopy for a z=0.45 absorption line system 



4.2. Comparative chemistry 

It is interesting to compare the molecular abundances in 
the z=0.89 galaxy with those measured for various sources 
in the Local Universe, to characterize the type of clouds 
and chemistry. For this purpose, we selected a sample of 
archetype sources with a large number of detected molec- 
ular species, usually covered by large spectral surveys. 
Molecular abundances in circumstellar envelopes around 
evolved stars, such as IRC±10216, were not included be- 
cause of the distinctive chemical segregation within the en- 
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velope. We preferred to use data obtained from a single 
team/telescope, for the sake of homogeneity and to mini- 
mize beam effects. 

Despite the low density and poor shielding from the in- 
terstellar radiation field, a variety of molecules is already 
present in Galactic diffuse clouds (Ay < 1): a bout 15 dif- 
feren t species have been detected so far (see e.g. lLiszt et al.l 
2008), and the limit of their chemical complexity remains 
unclear. We mention here that column densities in diffuse 
clouds are most often determined through mm-wave ab- 
sorption toward extragalactic continuum sources, similar 
to the work presented in this paper. 

Translucent clouds (2<Ay<5) are expected to have a 
richer chemistry, with about 30 molecules detected so far, 
including highly unsaturated and very reactive hydrocar- 
bo n chains . We a dopt the molecular abundances compiled 
bv lTurnerl (12000H 

Cold dark clouds are the sites of low-mass star forma- 
tion, where an intricate gas-phase ion-molecule interstel- 
lar chemis try takes place. We adopt the abundances deter- 
mined bv lOhishi et all (jl992ft toward TMC-1, which is a 
prototype of dark clouds. 

Hot and dense cores, associated with massive-star for- 
mation harbor the most complex interstellar chemistry. 
There, the chemistry is characterized by the evaporation of 
dust grains, releasing large mole cules in the gas phase. We 
adopt the abundances derived bv lNummelin et all (|2000f ) in 
the region Sgr B2N toward the Galactic center as reference. 

Finally, we also include in our comparison three ex- 
tragalactic sources for which molecular abundances have 
been derived for a significantly large number of species: 
the nuclea r regions of the sta rbu rst galaxies NGC 253 and 
NGC4945 (|Martm et al.ll2006l and lWang et al.ll2004 respec- 
tively), as well as the star forming region N1 13 in the Large 
Magellanic Cloud (XMC. IWan7et~aUl200l . 

The comparison of abundances derived in FG0.89SW 
with those observed toward these various sources is illus- 
trated in FigJHl Molecular abundance ratios, normalized 
to FG0.89SW and detailed by species, are also given in 
FigM Several species in our survey, such as HCO + , HCN, 
C2H, and C-C3H2 are commonly observed in Galactic dif- 
fuse clouds. More complex molecules (e.g. carbon chains) 
are only seen toward cold dark clouds or in warm and dense 
clouds. 

Comparing molecular abundances relies on the validity 
of various hypotheses. First, molecules have to be co-spatial 
and share the same source covering factor. The size of the 
SW and NE compact continuum components corresponds 
to a few pc in the plane of the z=0.89 galaxy. We assume 
a common source covering factor of unity for all molecules. 
However, the absorbing gas could be the mix of (most likely 
several) dense cores embedded in a more diffuse compo- 
nent within the continuum beam illumination. Our derived 
column densities and molecular abundances are thus lower 
limits. In addition, opacity effects, if not correctly taken 
into account, could easily lead to erroneous molecular abun- 
dances. For this reason, we d o not use the abundance s of 
HCO+, HCN, and C 2 H from iNummelin et al.l (|2000l ) to- 
ward Sgr B2N, as they are probably largely underestimated. 
Finally, a reference species needs to be assigned for normal- 
ization, in order to compare the relative molecular abun- 
dances. We choose to normalize the column densities to that 
of molecular hydrogen, as commonly done in the literature. 
For this, we assume a value N(H2)=2 x 10 22 cm~ 2 toward 



the SW component, as estimated by IWiklind fc Combesl 
(|l998l ). This value is consistent with the measurement of 
A(H)= (3.5 ± 0.6) x 10 22 cm" 2 obtained from X-ray absorp- 
tion bv lMathur fc Nairl (11997ft assuming that all absorption 
is due to the SW component. The absorb ing gas is indeed 
most ly molecular along this line of sight (|Chengalur et al.l 
1999)', so that A(H 2 ) should be close to 2x7V(H). The 
NE line of sight, on the other hand, intercepts the disk of 
the z=0.89 galaxy at a galactocentric distance of ~4 kpc. 
There, the extinction Ay is roughly one order of magni- 
tude lower than toward the SW image of the quasar (e.g. 
iFalco et alll 19991 see also FigfT]), and the column density of 
molecular gas is likely to be similarly lower. We assume a 
H 2 column density of 1 x 10 21 cm~ 2 toward FG0.89NE. This 
choice makes the fractional abundances similar to those in 
Galactic diffuse clouds (Tableland Fig|H|). While the nor- 
malization to H2 introduces an uncertainty in the absolute 
molecular abundances, abundance ratios for a given source 
are not affected, and should therefore bear greater signif- 
icance (again provided that the species are co-spatial and 
have an identical source covering factor). 

To compare the ov erall trend between the various 
sources, and inspired bv iMartrn et all (|2006t ). we have cal- 
culated an abundance estimator Xij > defined as the average 
of the logarithm of abundance ratios between two sources 
i and j 

x ^^g logl0 (iS)' (5) 

where JQ /j stands for the abundance relative to H 2 of 
molecule k in source i, and the sum runs over all N mo i 
molecules detected in common for the two sources i and 
j. A dispersion in the estimator was also calculated to re- 
flect potential large variations in the abundance ratios be- 
tween species. We do not take into account non-detections 
in this estimator. Abundance estimators were also calcu- 
lated in the same manner for the following sub-samples of 
molecules: carbon chains (i.e. with more than two chained 
carbon atoms, such as I-C3H, C4H, C-C3H2, I-C3H2, and 
HC 3 N), sulfur-bearing species (i.e. SO, NS, H 2 CS, C 2 S, and 
SO+), and saturated molecules (CH 3 OH and CH 3 NH 2 ), 
presumably all formed on dust grains. The value of Xij 
was calculated between FG0.89SW and the other sources 
only if the number of common species for two objects was 
>2. The abundance estimators are given in Table [TU] and 
illustrated in FigfTOl 

On the basis of the abundance estimators including all 
species, we obtain a basic sequence of increasing molecular 
abundances such as LMC < diffuse clouds < FG0.89NE < 
NGC253 < FG0.89SW < NGC4945 < TMC 1 < translu- 
cent clouds < Sgr B2N, indicating, unsurprisingly, that the 
chemical abundances in FG0.89SW are in-between those in 
typical Galactic diffuse and translucent clouds. 

The dispersions for Sgr B2N are large, though it is 
clear that sulfur-bearing and especially saturated species 
have much higher abundances in this molecular-rich ob- 
ject. In TMC 1 and translucent clouds, the abundance of 
carbon chains is globally a factor of a few higher than in 
FG0.89SW, while that of sulfur-bearing species is one order 
of magnitude higher. 

The overall similarity between the molecular abun- 
dances of FG0.89SW and the nuclea r region of the star- 
burst galaxy NGC253 is surprising. IMartm et all (|2006h 
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found good agreement between the molecular abundances 
in NGC253 and those at the position Sgr B2(OH) in the 
Galactic center molecular cloud Sgr B2. Shocks are thought 
to play an important role in the heating and chemistry to- 
ward these nuclear clouds, triggering the disruption of dust 
grains, and releasing the products of grain-surface chem- 
istry into the gas phase. It is tempting to invoke simi- 
lar chemical processes in FG0.89SW, as they could ex- 
plain all of the similarities between the molecular abun- 
dances, the high kinetic temperature of the ga s, and the rel- 
ative l y large linewidth (s ee also discussions in lHenkel et al.l 
l2008tlMenten et al.|[20Q8t ). Moreover, the elemental compo- 
sition of the gas in both galaxies is expected to be dom- 
inated by the nucleosynthesis products of massive stars, 
because of their shorter timescales. The major difference, 
though, is that FG0.89SW clouds are located in a spiral 
arm at a galactocentric distance of about 2 kpc, and gas 
densities are likely lower. Abundances in NGC4945 show 
some differences, albeit small, when compared to those in 
NGC253. They are interpreted as different starburst evolu- 
tionary states ()Wang et all l2004t [Martin et alj|2006l) . The 
systematically lower molecular abundances in the LMC- 
N113 region are attributed to a photon domi nated region 
(PDR ) in a nitrogen deficient environment (jWang et al.l 
2009). Nevertheless, we emphasize that, because of the lim- 
ited angular resolution of these extragalactic observations, 
the corresponding molecular abundances reflect the average 
of the different types of cloud gas properties and excitation 
conditions within the beam. This is much less the case for 
absorption line studies toward quasars, where the spatial 
resolution is set by the small size of the background con- 
tinuum emission. 



4.3. Kinematics 

4.3.1. Velocity centroids and line profiles 

The limited velocity resolution of our survey (6-10 kms -1 , 
although slightly improved after the Doppler correction) 
prevents us from conducting a detailed analysis of the kine- 
matics and line profiles. Nevertheless, the high S/N of most 
of the lines allows us to determine the velocity centroids of 
a large number of species from Gaussian fitting, with an un- 
certainty of the order of one kms -1 or better. Inter-species 
velocity offsets can then be tested statistically. The results 
are shown in FigfTTl where we show the FWHM versus the 
velocity centroid derived from Gaussian fits for species seen 
toward the SW absorption. 

Excluding the strong and likely saturated absorption 
of HCO + and HCN, we obtain an average velocity of 
-1.78 ± 1.83 kms" 1 toward FG0.89SW. Few species, 
namely CH 3 NH 2 , CH 3 OH, 1-C 3 H, HNCO, SO+, and HC 3 N, 
have velocities that differ from this average by more than 
1 (T. These kinematical differences might be due, e.g., 
to chemical segregation along the line of sight, different 
source covering factors per species, excitation conditions, 
asymmetric line profiles, or a limited velocity resolution 
and/or S/N of the observations. Another interpretation is 
explored in the next section. After removing the outliers 
listed above, we obtain a new estimate of the average ve- 
locity: — 1.51 ± 0.73 kms -1 , the dispersion dropping by a 
factor of more than two with respect to the first value in- 
cluding outliers. 



It will be interesting to investigate the chemical com- 
plexity toward FG0.89SW when searching for more species, 
especially light hydrides and larger molecules. The molec- 
ular absorber located at z=0.68 in front of the quasar 
B 0218+357 would also be a suitable target for these in- 
vestigations. 



o Carbon chains 

S— bearing species 
□ Saturated species 



1 - 
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Fig. 10. Comparison of the abundance estimators Xij as 
defined in EqjSl for different sources and sub-samples of 
molecular species (see text), with respect to FG0.89SW. 



30 



25 



E 20 

CD 



15 




Velocity (km s ) 

Fig. 11. Diagram of the FWHM versus velocity centroids, 
derived from Gaussian fits to the lines toward FG0.89SW. 
The average velocity (full line) and 3<7 dispersion (dashed 
lines), calculated from 22 species (see text), are indicated. 
The two points with a linewidth larger than 25 kms -1 cor- 
respond to HCO + and HCN, the lines of which are likely 
saturated. 
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4.3.2. Constraints on the variation in the fundamental 
constants 

We now discuss the interpretation of the velocity offsets in 
terms of variations in the fundamental constants of physics, 
and more particularly, in the proton-to-electron mass ratio. 

Searches for observational evidence of any secular vari- 
ation in the fundamental constants (FC) of physics has at- 
tracted much attention lately. Any discovery of such evi- 
dence would indeed contradict the Standard Model in the 
invariance principle, and open new horizons in physics re- 
search. We note that several theories predict that the fun- 
damental constants vary, an d we r e fer th e interested reader 
to a review of this field bv lUzanl (|201ll ). A description of 
the various methods that have been applied to constrain 
the variation in the fundamental constants to date, their 
results, and potential systematics are also given in this re- 
view. 

Spectroscopic observations of absorption lines toward 
high-redshift quasars naturally offer very long time lapse, 
possibly longer e.g. than the age of the Solar System, that 
could allow even small variation rates to be measured. 
While different methods have been developed, they all con- 
sist of a change in the frequency of atomic or molecu- 
lar transitions induced by a potential variation of the (di- 
mensionless) fine structure constant (a), proton-to-electron 
mass ratio (/i=m p /m e ), and nuclear g-factor (g n ). One can 
write this change as 



v a (i 



Ag n K AC lc\ 

K 9— = 2^ k c-tt, (6) 



where the Ki coefficients reflect the sensitivity of a change 
in the frequency of the transition v to a change in the con- 
stant C=a, [i and/or <?„. The comparison of two different 
lines is then given by 



AV 

c 



V, - V, 



E 



(K Cj - K c ,i) 



AC 

C ' 



(7) 



where AV^ is the velocity difference between the lines i 
and j induced by changes in FCs. From EqfTl it is clear 
that only the comparison of lines with different Kq sensi- 
tivity coefficients yields a possibly non-null velocity shift. 
For example, one has |iC a |<<l, K^—— 1 and |iC 5 |<<l for 
rotational transitions. For hyperfme transitions, K a =— 2, 
K^=— 1, and K g =— 1. Therefore, the comparison of hyper- 
fine and rotational lines is sensitive to a variation of AF/F, 
where F=a 2 g n . 

Inversion lines of ammonia have a high sensitivity to any 
variation in /u, with if,, =—4.4 6, but are nearly insensitiv e 
to changes in other constants (|Flambaum &: Kozlov 2 007f ). 
As a result, the comparison of NH3 inversion lines and other 
rotational transitions provide a test of ^ . This method was 
used toward the quasars B 0218+357 and PKS 1830-211 
(the only two known high-z ammonia absorbers) and yields 
no significant evidence of changes in fi, with constraints on 
^ d own to a few 10" 7 (jMurphv et alj|2008t IHenkel et all 
2009; iKanekarl l201lh , corresponding to a velocity offset 
of the order of a few tenths of kms -1 . Nevertheless, the 
method is limited by the following systematic effects: (1) 
local velocity offsets between species; (2) non-LTE and sat- 
uration effects for NH3 hyperfine strucure (if hyperfme com- 
ponents are unresolved or poorly detected) and optically 



thick lines, respectively; (3) time variations (if observations 
are taken at different epochs); and (4) changes in the con- 
tin uum morphol o gy at significantly different frequencies. 

IHenkel et al.l (|2009[ ) measured a velocity of —3.70 ± 
0.43 kms -1 for HC3N, consistent with our value of — 3.90± 
0.46 kms -1 . They used this value as a reference to constrain 
the velocity shift of NH3 inversion transitions, obtaining a 
3 a upper limit on of 1.4xl0~ 6 . Taking our average 
velocity (without outliers) as reference instead, we derive 
a new value ^ = (—2.04 ± 0.74) xl0~ 6 from the ammonia 
velocity shift, assuming no time variations between the dif- 
ferent observations and that the shift is indeed due to vari- 
ations in ji. Taken at face value, the velocity dispersion 
of the different species in our survey introduces a signifi- 
cant contribution as kinematical noise, potentially limiting 
constraints on — to ~2xl0~ 6 using inversion lines of am- 
monia. 

Interestingly, methanol was also found to have tran- 
sitions with high sen sitivity to change in the proton-to- 
electron mass ratio (pansen et al.l 1201 lh . The crude off- 
set between the velocity of methanol and the average ve- 
locity is AV=— 3.75 ± 0.89 kms -1 . Taking the sensitiv- 
ity coefficient K u =— 7.4 for the m ethanol transition at 
60.531477 GHz (jJansen et al.ll201lD and provided that the 
velocity offset is due only to a variation in /i, we obtain 
^ = (-1.95 ± 0.47) xl0~ 6 , that is, a statistical 4 a sig- 
nificant detection of a variation in /i at z=0.8858fl 111 
agreement with the ammonia result. However, since we de- 
tect only one methanol line in our survey, we consider this 
result as only tentative. More conservatively, our obser- 
vations give an upper limit |-^|<4xl0~ 6 , which, impor- 
tantly, takes into account the velocity dispersion of a large 
number of molecular species. More species and lines with 
different sensitivity coefficients should be observed at high 
velocity resolution and high S/N to check th e comp atibility 
with variations in /1. We note that IKanekarl (|201 lh derived 
^ = (-3.5 ± 1.2) xlO -7 at z=0.68 from the comparison of 
NH3 inversion lines and rotational lines of two species, CS 
and H 2 CO, seen in absorption toward B 0218+357. 



5. Summary and conclusions 

We have completed the first unbiased spectral survey 
toward an intermediate-redshift molecular absorber, 
the z=0.89 galaxy located in front of the quasar 
PKS 1830—211. Observations, performed with the 
Australian Telescope Compact Array, cover the 7 mm 
band from 30 to 50 GHz (equivalent to the rest-frame 
frequency interval 57-94 GHz), and reach a sensitivity of 
the order of a few times 10~ 3 of the continuum level, at a 
spectral resolution of 1 MHz. 

The brightness of the background quasar, the large col- 
umn density of absorbing molecular gas, the large instan- 
taneous instrumental bandwidth, the wide frequency cov- 
erage of our survey, and the frequency band most sensitive 
to low energy transitions of numerous molecules, all concur 
in allowing us to achieve the first large chemical inventory 
in the interstellar medium of a distant galaxy. 

From these observations, we derive the following results: 



5 That is, a look-back time of 7.24 Gyr for 
#o=71 kms" 1 Mpc" 1 , n A f=0.27, and fi A =0.73. 
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— We detect a total of 28 molecular species toward the 
SW image of the quasar (of which more than half are 
detected for the first time toward this target and eight 
are detected for the first time outside the Milky Way: 
SO+, I-C3H, I-C3H2, H 2 CCN, H 2 CCO, C 4 H, CH 3 NH 2 , 
and CH3CHO), and six toward the NE image. 

— The rotation temperatures are close to the expected 
cosmic microwave background temperature, 5.14 K, at 
0=0.89. This absorber, with strong absorption along two 
independent lines of sight, offers a unique opportunity to 
establish a robust and accurate measurement of Tqmb 
at z=0.89. 

— The fractional molecular abundances toward the SW 
absorption are found to be in-between those in typical 
Galactic diffuse and translucent clouds. 

— Toward the NE absorption, where the column density 
of molecular gas is about one order of magnitude lower, 
the fractional abundances are comparable with those in 
diffuse clouds. 

— The isotopic ratios of C, N, O, and Si at 7 mm confirm 
and comple ment the C, N , O, a nd S ratios previously 
obtained by I Muller et al.l (|2006l ) at 3 mm. Significant 
differences are found from the solar and local ISM val- 
ues, which could be explained by the limited timescale 
for nucleosynthetic processing at z=0.89. Deuterated 
species remain undetected down to a D/H ratio of 
<7xl0- 4 . 

— We discover several new, but weak, velocity compo- 
nents, at V=— 60, —224, and —300 kms" 1 , seen in the 
HCO+ and HCN (1-0) lines. 

— The analysis of the velocity centroids yields a disper- 
sion of nearly 2 kms -1 toward FG0.89SW. Among 
the species showing a significa nt offset from the aver- 
age velocity, NH 3 (observed bv lHenkel et al.l l2009) and 
CH3OH both have a high sensitivity to variations in 
the proton-to-clcctron mass ratio /i. We derive a con- 
servative upper limit of |-^|<4xl0~ 6 at z=0.89, which 
includes the velocity dispersion (hence the uncertainty 
caused by potentially different spatial distributions) of 
a large number of molecular species. Further observa- 
tions of lines with different sensitivity coefficients will 
be required to place stronger constraints on variations 
in [i. 

The new generation of instruments currently deployed 
(ALMA, EVLA, NOEMA) will offer new capabilities in 
terms of frequency coverage, instantaneous bandwidth, sen- 
sitivity, and spectral resolution that will be of great benefit 
to molecular absorption studies. It will be interesting to in- 
vestigate the chemistry and its complexity in this distant 
z=0.89 galaxy, as well as toward other known molecular 
absorbers. 
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Fig. 2. a) Full 30-50 GHz spectral band observed toward PKS 1830—211 with the ATCA, normalized to the continuum 
flux of the SW image of the quasar (38% of the total flux), b) Same spectral interval, but limited to a scale of a few 
percent of the SW continuum intensity, to illustrate the noise level over the band. The rest-frame frequency axis is given 
on top of the figure. 



Table 3. Census of species detected toward FG0.89SW. 
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Notes. In parenthesis, species detected in other studies; bold face, new detections; underlined, first extragalactic detections. 

O including H 13 CO + , HC ls O + , and HC 17 + ; A including H 13 CN and HC 15 N; including HN 13 C and H 15 NC; t including 29 SiO. 



Re ferences. (a)IWiklind fc Combed (Il996bfl: (b ) lGerin et all (|1997l ): fcl lWiklind fc Combesl (Il998f) : (d ) IChengalur et all 1119991) : 
(e) iMenten et all (119991 ): (p iMuller et al.1 (I2006D . including isotopologues: (s) iHenkel et al.1 (|200W : (h) iMenten et al.1 l|2008f)~ i') 
iBottinelli et alJ l|2009l ): ffl lHenkel et all l|2009fl 
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Table 4. Molecular data for species detected toward the SW absorption. 
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Notes. (a) For species with ortho/para symmetry, the two forms were considered as distinct species and the partition function 
Q was thus calculated accordingly, with the weight factors g=3 for ortho form, and g=l for para form. Hyperfine structure was 
taken into account when necessary. 

' 6 - ) Number of lines detected or resolved (in the case of hyperfine structure). 

' c ' Rotation temperature and column density obtained from Monte Carlo simulations of rotation diagrams. Uncertainties 
correspond to a 95% confidence interval. 

Column density, assuming a rotation temperature of 5.14 K and LTE. 
^ Possibly underestimated due to the large opacity of the line. 
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Table 5. Line parameters and fitting results for the NE absorption component. 



Species 


Vb 


AV 
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H 13 CN 
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38.1 (0.0) 
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Notes. Negative values (except for velocities) represent either the (3 a) upper (for Nlte) or lower (for [SW]/[NE]) limits. 
Table 6. Upper limits (3a) for non-detections of some other interstellar species toward FG0.89SW. 
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■6 
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■8 
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/ 
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Fig. 8. Comparison o f relative molecular abundances in FG0.89SW (from most to less abundant), diffuse cloud s 
(iLucas fc Lisztll2000al|bl: iLiszt fc Lucasl 120011: iLucas fc Lisztl l2002t ITiszt et al.ll2004h. translucen t clouds ([Turner! [2000h . 
TMC-1 (lOhishi et al.llT992ft . Sgr B 2(n1 (iNummelin et alJl2000D , LMC -N113 (iWang et aUl2009lV and the nuclear region 
of the starburst galaxies NGC253 (jMartm et al.ll2006t ) and NGC4945 (|Wang et alJl2004T >. 
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Table 9. Relative abundances [X]/[H 2 ] (10~ n ). 



Species 


FG0.89SW 


FG0.89NE 


Diffuse 
clouds 


Translucent 
clouds 


TMC 1 


Sgr B2(N) 


LMC 
N113 


NGC253 


NGC4945 


C 2 H 


6242.4 


2924.5 


2600.0 


6600.0 


7500.0 


- 


- 


2000.0 


12600.0 


HCN 


1520.4 


733.0 


300.0 


3600.0 


2000.0 


- 


30.0 


500.0 


2000.0 


H 2 CO 


1196.3 


422.4 


400.0 


630.0 


2000.0 


- 


80.0 


250.0 


200.0 


HCO+ 


875.4 


746.7 


200.0 


200.0 


800.0 


- 


20.0 


160.0 


- 


CH3OH 


849.9 


- 


- 


1800.0 


200.0 


20000.0 


50.0 


1260.0 


790.0 


C 4 H 


656.7 


- 


- 


1700.0 


2000.0 


- 


- 


- 


- 


HNC 


509.2 


116.8 


60.0 


250.0 


2000.0 


- 


10.0 


100.0 


3200.0 


CH3CCH 


300.0 


- 


- 


270.0 


600.0 


400.0 


- 


630.0 


400.0 


C-C3H2 


264.5 


125.0 


140.0 


3600.0 


1000.0 


3.0 


20.0 


50.0 


500.0 


CH 3 NH 2 


160.6 


- 


- 


- 


- 


30000.0 


- 


- 


- 


SO 


128.1 


- 


150.0 


3200.0 


500.0 


2000.0 


80.0 


130.0 


100.0 


N 2 H+ 


116.6 


-33.4 


-0.4 


100.0 


50.0 


- 


3.0 


- 


60.0 


H2CCO 


88.3 


- 


- 


110.0 


100.0 


20.0 


- 


- 


- 


CH 2 NH 


74.5 


- 


- 


1500.0 


- 


1000.0 


- 


80.0 


- 


HCO 


65.0 


- 


- 


- 


- 


- 


- 


280.0 


- 


CH 3 CN 


65.0 


- 


- 


-110.0 


- 


3000.0 


- 


30.0 


- 


CH3CHO 


62.0 


- 


- 


1100.0 


60.0 


20.0 


- 


- 


- 


NS 


57.5 


- 


- 


- 


- 


1000.0 


- 


60.0 


- 


HC 3 N 


55.9 


- 


- 


50.0 


600.0 


500.0 


3.0 


60.0 


100.0 


H 2 CS 


37.6 






240.0 


300.0 


2000.0 




60.0 




SiO 


37.5 




10.0 


10.0 


-0.3 






13.0 




I-C3H 


35.7 






500.0 


50.0 










HNCO 


31.1 






240.0 


20.0 


60.0 




160.0 


250.0 


C 2 S 


24.2 






240.0 


800.0 






20.0 




HOC+ 


15.9 




5.0 










30.0 




H 2 CCN 


14.6 








500.0 


3.0 








I-C3H2 


10.5 






200.0 












SO+ 


8.4 






100.0 













Notes. Negative values represent an upper limit. H2 column densities of 2xl0 22 and lxlO 21 cm 2 were assumed for FG0.89SW 
and FG0.89NE, respectively. 

Refere nces. D i ffuse clouds: iLucas fc Lisztl (l2000al|H) ; iLiszt fc Lucas! (12001T); ILucas fc Lisztl (|2002f ) ; ILiszt et all (I2004D ; Tran slucent 
clouds: iTurneil (|2000f ): TMC-1: lOhishi et all dl992Tl: : SgrB2CN): iNummelin et al.1 (120001 ): LMC-N113: IWang et al.l <|2009l 1. using 
N(R 2 )=2 x 10 22 cm -2 ; NGC2537 lMartm et al.l lj2006f l: NGC4945 TWang et al.l l|2004h . 

Table 10. Abundance estimator. 



Source 


All 


Carbon 




S-bearin 




Saturated 




species 


chains 




species 




species 


FG0.89SW 


















LMC N113 


-1.3( 0.4) [ 9] 


-1.2(0.1) 


2] 


-[1] 




-[1] 


Diffuse clouds 


-0.5( 0.3) [ 9] 


-[1] 




-[1] 




-[0] 


FG0.89NE 


-0.4( 0.2) [ 6] 


-[1] 




-[0] 




-[0] 


NGC253 


-0.1( 0.4) [19] 


-0.3( 0.4) 


2] 


0.0( 0.1) [ 


4] 


-[1] 


NGC4945 


0.1( 0.4) [11] 


0.3( 0.0) [ 


2] 


'-[1] 




-[1] 


Transl. clouds 


0.5( 0.7) [22] 


0.8( 0.5) [ 


5] 


1.1( 0.2) [ 


4] 


-[1] 


TMC1 


0.4( 0.6) [19] 


0.6( 0.3) [ 


4] 


1.0( 0.4) [ 


3] 


-[1] 


Sgr B2N 


0.5( 1.1) [13] 


-0.5( 1.4) 


2] 


1.4( 0.2) [ 


3] 


1.8( 0.4) [ 2] 



Notes. The abundance estimators were calculated according to Eqf5] The dispersion is given in parenthesis and the number of 
molecules detected in common between FG0.89SW and other sources for each sample is given within brackets. 
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Fig. 9. Ratios of molecular abundances of FG0.89SW to other sources, detailed by species. 
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Appendix A: Spectral survey and line catalog 

The full spectrum, from 30 to 50 GHz, is available on-line as 
an ASCII file at the Strasbourg Astronomical Data Center 
(CDS). 
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Fig. A.l. Spectrum of the ATCA 7 mm survey toward PKS 1830—211. The grey envelope indicates the 3 a noise level. 
The LTE synthetic model, assuming T rot =b.lA K for all molecules, is overlaid in red. The spectrum in magenta, offset at 
a level of 0.995, gives the relative differences between observations taken in September 2009 and March 2010, normalized 
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Table A.l. Parameters and fitting results for lines detected toward 
FG0.89SW. 
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Table A.l. Continued. 



Rest Freq. 


Transition 
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Notes. For species with hyperfine structure, the relative intensities were fixed, and the given uncertainty corresponds to that in 
the total integrated intensity. 
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Appendix B: Notes on individual molecules 

In this section, we give short reviews of the different species 
and discuss constraints or new results obtained from this 
survey. 

B.l. C 2 H 

The highly reactive ethynyl radical was fir st discovered 
and id entified toward the Orion nebula by I Tucker et al.l 
(|1974[) based on its hyperfine structure, while it had still 
not been observed in the laboratory. C2H has turned out 
to be both an ubiquitous and abundant speci es in the inter- 
stellar medium, including in diffuse clouds (jLucas fc Lisztl 
2000b]) where its abundance is of the order of a few 10 -8 
with respect to H 2 . 

The ground state transition N=l-0 at ^0—87 GHz is 
decomposed into six components by hyperfine splitting. 
Toward PKS 1830—211, the strong est hyperfine c ompo - 
nents were previously observed by iMenten et al.l (l999) 
with the VLA, although at low spectral resolution. We de- 
tect here all six iV=l-0 hyperfine components arising from 
the SW absorption, as well as the two strongest hyperfine 
components (the J=3/2-l/2, F=2-l and -F=l-0) related to 
the NE absorption at -147 kins" 1 (see Fig fB.ll) . The ab- 
sorption depths reflect the LTE relative intensities of the 
hyperfine components, indicating that the lines are opti- 
cally thin. 




46.3 46.35 46.4 

Sky Frequency (GHz) 



Fig. B.l. Spectrum of the C 2 H N=l-0 line. The different 
hyperfine components are indicated with their relative in- 
tensities for the SW (top) and NE (bottom) absorptions. 
Note the clear detection of the two strongest hyperfine 
components arising from FG0.89NE. The red curve shows 
the LTE synthetic spectrum, taking column densities of 
1.3xl0 15 and 4xl0 13 cm" 2 toward the SW and NE im- 
ages respectively, and a rotation temperature of 5.14 K. 



two lines are likely saturated and their opacities could be 
underestimated. In particular two findings seem to indicate 
that the HCO + column density might be underestimated 
by a factor of 2-3. First, the HCO+/H 13 CO+ ratio toward 
FG0.89SW is a factor of ^2 lower than the corresponding 
HCN/H 13 CN and HNC/HN 13 C ratios. Second, the ratio 
of the column densities along the two lines of sight, 
SW/NE is also a factor 2-3 lower for HCO + than for other 
species (see Table [5]). In any case, the uncertainty in the 
opacity introduced by the uncertainty in the continuum 
background illumination is nevertheless large for both 
HCO+ and HCN (see §Q 

Time variations and the discovery of additional velocity 
components are discussed in § 13.21 and § 13. 51 respectively. 
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Fig.B.2. Full spectra of the HCO+, HCN, and HNC 
1-0 transitions, showing the absorption toward the SW 
(0 kms -1 ) and NE (—147 kms -1 ) images of the quasar. 
A zoom in with a narrower intensity scale, shown in Fig(SJ 
reveals additional velocity components. 



B.2. HCN, HNC and HCO+ 

Molecular absorptio n was first rep o rted toward 
PKS 1830-211 by iWiklind fc Combes! (11996bft after 
the detection of the HCN, HNC, and HCO+ 2-1 transi- 
tions, allowing the redshift of the intervening galaxy to be 
derived. FigEHshows the full spectra of the HCO+, HCN, 
and HNC 1-0 transitions observed in our survey. The first 



B.3. N 2 H+ 

N 2 H + is thought to be one of the most abundant interstel- 
lar ions, and plays a key role in nitrogen chemistry. It is 
detected in a large number of G alactic clouds. Howe ver, it 
remains elusive in diffuse clouds (jLiszt fc Lucasl[200lT) , with 
an abundance upper limit more than two orders of magni- 
tude lower than its abundance in FG0.89S W. We note that 
N 2 H + is also rare in translucent clouds, e.g.. lTurnerl (|1995al) 
having detected it in only 2 out of 16 such clouds. 
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B.4. HOC+ 

The hydroxymethylidynium ion, HOC + , isomer of the 
formyl ion, was fir st tentatively identified in Sgr B2 by 
IWoods et al.l ([19831 ). Its presence in several dense molec- 
ular clouds was then c o nfirme d bv lZiurvs fc A pponi (jl995h 
and lApponi fc Ziurvsl ([19971 ). with a relative abundance 
HCO+/HOC+-360-6000. This abundance ratio is how- 
ever about one order of magnit ude lower (—70-120) in 
diffuse clouds (jLiszt et al.l I2004D. and se e mingly has an 
in-between value in PDR: lApponi et all (119991) derived 
HCO+ /HOC+-270 toward the Orion bar iMartm et al.l 
(2009) find a ratio -70 toward NGC253, similar to that 
in diffuse clouds. 

In FG0.89SW, we derive a ratio HCO+/HOC+ of -60, 
although this value might be slightly underestimated be- 
cause of the large opacity of HCO+. The H 13 CO+/HOC+ 
ratio is indeed 2.6±0.1. Assuming a 12 C/ 13 C ratio of 35, the 
HCO+/HOC+ ratio in FG0.89SW is comparable to that in 
diffuse clouds. 



Turner et al.lll999f). but is not det ected in cold dark clouds 
Turner! 1 199H iDickens et al.lll997fl. It ha s been tentatively 



detected in NGC253 (IMartm et alll2006D and has been ob- 
served in Arp 220 (jSalter et aLll2008l) . Methanimine is a 
prebiotic molecule, precursor to more complex molecules 
such as glycine. 

B.7. CH 3 NH 2 

Methylamine is the saturated and terminal product of the 
hydrogenation series based on the cyanide radical: CN- 
HCN-CH 2 NH-CH 3 NH 2 . 

Methylamine has a very complex spectrum resulting 
from the combination of internal rotation of the CH3 group, 
the asy mmetry of the molecule and the inversion of the N H 2 
group ([Takagi fc Koiimal 119731: lllvushin fc Lovasll2007l) . It 



was first det ected in Sgr B2 bv lKaifu et al. ( 19741 ) (see also 
lTurnerlll99ll ). No extragalactic detection has been reported 
to date. 



B.5. HCO 

The formyl radical, HCO, was d i scover ed in the inter- 
stellar medium by ISnvder et al.l ([1976D and has been 
identified as a tracer of the clo u d-illuminated inter - 
face in PDRs (Ide Jong et all 119801: ISchilke et all l200lh . 
iGarria-Burillo et al.l (l200lT revealed widespread HCO 
emission in the nuclear starburst of M82, with HCO, CO, 
and HII emissi on in nested rings, HCO extending farther 
out in the disk. iGerin et all i|20()9i demonstrated that the 
HCO emission in the Horsehead nebula traces the PDR 
front, coincident with PAH and hydrocarbon emission. The 
abundance of HCO in the PDR is comparable to that of 
HCO + , whereas it becomes about one order of magnitude 
less in the dense regions shielded from UV radiation. To 
the best of our knowledge, HCO has not been detected in 
diffuse or translucent clouds to date, although it is located 
close in frequency to the H 13 CO+ (1-0) and SiO (2-1) lines. 

We clearly detect three (out of the four) hyperfine com- 
ponents of the ground state transition loi — 0oo 7 the re- 
maining low intensity component being buried in the noise 
level. The relative optical depths are consistent with LTE 
intensities, suggesting that the lines are optically thin. We 
derive an abundance ratio HCO + /HCO of —13, similar to 
the values observed in dense shielded regions. This value is 
even larger if we consider that the HCO 4 " opacity might be 
underestimated by a factor 2-3. 



B.8. CH 3 OH 

Interstellar methanol was fir st detected toward Sgr B2 and 
Sgr A by iBall et all ([1970D and is widely observed to- 
ward st ar forming re gions, dark clouds, and also translucent 
clouds ((Turner 1998). I t has been detect ed in various extra- 
galactic sources fe.g. lHenkel et al]|1987l ). Methanol is about 
one order of magnitude less abundant in cold/translucent 
clouds than in star forming regions, where its abundance 
can reach a few 10~ 8 relative to H 2 . It is not detected in 
diffuse clouds, down to upp er limits of a few lower than its 
abundance in cold clouds ([Liszt et al.l [2008). Methanol is 
the terminal product of hydrogenation of CO. 

The methanol spectrum has a multitude of transitions, 
of which only one is strong enough to be detected within our 
frequency coverage. Therefore, no excitation analysis can be 
deduced. Assuming LTE equilibrium with the CMB radia- 
tion field, we derive a relatively high methanol abundance 
of a few 10~ 8 relative to H 2 in FG0.89SW. Several strong 
lines (with predicted r- 0.02-0.4) near rest frequency of 300 
GHz and redshifted in the 2 mm band are (simultaneously) 
observable (e.g. with ALMA or PdBI) and could provide a 
stringent measurement of the methanol rotation tempera- 
ture. 

The calculations of l Jansen~e t al. (201l|) have shown that 
methanol could be a sensitive probe of spatial and temporal 
variations in the proton-to-electron mass ratio. 



B.6. CH 2 NH 

We identify the absorption feature near 33.93 GHz with the 
loi — 0oo ground state transition of methanimine, as it is the 
simplest interstellar molecule with lowest energy transition 
close to that frequency. The hyperfine structure is not re- 
solved. Only one line of methanimine is located within the 
limits of our frequency coverage. The strongest transitions, 
with predicted r > 0.1 based on the loi — 0oo absorption 
depth, have rest frequencies of 127.9, 166.9, 225.6 and 284.3 
GHz. Of these, only the second and last fall in an ALMA 
band for a redshift of z=0.89. 

Methanimi n e wa s first detected in Sgr B2 by 
iGodfrev et all (|1973l ) and has been further observed 
in several star-forming regions and translucent clouds 



B.9. I-C3H 

We identify the feature around 40.4 GHz with the J=7/2- 
5/2 transition of the 2 Ii\/2 rotational ladder of the propy- 
nylidyne radical (^o=76.2 GHz). The transition is split 
into two A-doublets (see Fig lB.3[) . approximately equal in 
strength and spaced by about 5 MHz in the rest frame, 
which secures our identification. The hyperfine structure is 
unresolved. No other transitions of I-C3H are detected in 
our frequency coverage. 

Linear propynylid yne was first identified by 
iThaddeus etail (Il985al) toward IRC+ 10216 and TMC 1 . 
It is observed in translucent clouds dTurner et al. 2000) 
but re mains undetected in diffuse clouds (|Lucas fc Lisztl 
l2000bl although this might just be a sensitivity issue, as 
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their upper limits of 1-C3H/HCO+ is still higher than the 
value we derive). No extragalactic detections have been 
reported to date that we are aware of. 

The cyclic isomer form C-C3H is d etected in dense 
( Yamamoto et al. 1987t iFosse et all 1200 1|) and translucent 
(jlurner et al.ll2000l) clouds, with C-C3H/I-C3H ratio ~5. No 
transitions of C-C3H are detected in our survey. In particu- 
lar, the non-detection of the N KaKc=^n-^\\ transition at 
a rest-frame frequency of 91.5 GHz, which is connected to 
the ground state, yields an upper limit of ~3xl0 12 cm~ 2 , 
and a C-C3H/I-C3H ratio <0. 4. We note that a lo w ratio 
c-C 3 H/l-C 3 H< 1 was found bv lTurner et~aTl (|2000D toward 
one translucent cloud (CB228). 



opacity of r=0.007, 0.013, and 0.006, respectively) could 
confirm the presence of HCS + . 




40.39 40.4 40.41 40.42 

Sky Frequency (GHz) 

Fig. B.3. Spectrum of the 1-C 3 H J=7/2-5/2 transition. 
The locations of the two A-doublets with hyperfine struc- 
ture are indicated. 



B.10. C-C3H2 

The ring molecule cyclo-propenylidene C- C3H2 is a fairly 



ubiqui tu ous 
1985b!: 



interstellar 



Thaddeus et all 
~ 1981 



spec ies ( e.g., 

.Matthews k Irvine! 119851: ICox et al 
Madden et al.lll989l ). 

We detect 2 ortho and 2 para lines, yielding an or- 
tho/para ratio of 2.8 ±0.1, comparable to the statistical 
weight 3:1. 

A tight correlation exists betwe en the abundances of 
C2H and C-C3H2 in diffuse clouds (jLucas k Liszt! [200011 
lOerin et al.ll20ll . with A(C 2 H) = (28±1.4)iV(c-C3H 2 ). We 
derive a comparable abundance ratio C2H/c-C3H2=23.6 ± 
0.2. 

The excess absorption in the blue wing of the C-C3H2 
2i 2 -loi line (see Fie lB.4p could be due to the known veloc- 
ity compo nent at ~-20 km s" 1 toward FG0.89SW (see e.g. 
Fig. 5 from lMuller et aill2006f) . Alternatively, we note that 
this feature could be due the HCS + 2-1 line, which would 
fall at ^—32 km s -1 with respect to the center of the C-C3H2 
2i 2 -loi line. HCS + is indeed detected in diffu se clouds with 
an ab undance of ~3xl0~ 10 relative to H2 (jLucas k Liszt! 
2002). If its abundance were similar in FG0.89SW, then 
the LTE predicted absorption depth would be of the order 
of r=0.015, i.e. compatible with the observations. There 
are no other lines of HCS + within our frequency coverage. 
Observations of the J=l-0, 3-2, or 4-3 transitions (red- 
shifted to 22.6, 67.9, and 90.5 GHz, and with predicted 
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Fig. B.4. c-C3H 2 transitions. The dashed curve on the wing 
of the 2i2 — loi line corresponds to the LTE prediction 
for the HCS + J=2-l line with a HCS + column density of 
cm" 



B.ll. /-C3H2 

Propadienylidene, I-C3H2, is the linear isomer o f C-C 3H2 
that was first discovered by ICernicharo et all (|1991l) in 
the carbon-chain rich dark cloud TMC 1. Noticeably, 
iMaier et all (|201lh identified I-C3H2 as the probable car- 
rier of two conspicuous diffuse interstellar bands (at 4881 
and 5450 A). 

The ortho/para ratio is 2.4^q °, in agreement with the 
statistical weight of 3:1. 

The linear form is usually less abundant than the cyclic 
one, with a cyclic/linear ratio of 3-5 in diffuse clouds , 70 in 
TMC 1 and >150 in Sgr B 2 (|Cernicharo et al.lfl99ll fl999l 
see also lTurner et alll2000l ). We derive a cyclic/linear ratio 
of 25 ± 3, which is in-between the ratio observed toward 
diffuse clouds and TMC 1. 



B.12. H 2 CCN 

The non-linear cyanomethyl radical has a very complicated 
spectrum owing to the unpaired electron and hyperfine 
splitting induced by the H and N atoms. We detect two 
absorption features, corresponding to NK a ,K b =3o3-2o2 and 
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4o4-3o3- Both are from the ortho species (ortho states are 
ch aracterized by K a ev en and para states by K a odd). 

Ilrvine et ail (|1988l) first identified H 2 CCN toward 
Sgr B2 and TMC 1. It was a lso detected in the cir cumstel- 
lar envelope of IRC+10216 (jAgundez et al.ll2008l) . but, to 
the best of our knowledge, it has not been observed toward 
other sources. In particular, this is the first extragalactic 
detection of this molecule. 



B.13. CH 3 CCH 

Meth yl acetylene wa s first detected in Sgr B2 and Orion 
il (119731) . It is observed in dark clouds 
19811 ) and has been detected in the 
NGC4945 



in 

and M82 



A by Snvdcr & Bu 
(e.g., Ilrvine et al.l 
near by starburst galaxies NGC253 
(e.g. iMauersberger et alJll99ll: IWang et al.ll2004h . 

We detect the Jk—^q-^o transition, but not the Jr-=4 - 
3o, suggesting that the rotation temperature is higher than 
5.14 K. The kinetic temperature was fixed to 50 K (see 
lB~T4l) . 

We note that CH3CCH is the molecule in our survey 
which shows the miminum scatter in relative abundances 
between the different sources. 



Both are clearly detected. Remarkably, the if-ladder tran- 
sitions are detected up to K= 



=3 (see Fig lB.5p . with lower 
energy state up to 70 K above the ground level. We note 
that CH3CN is the only species in our survey showing ab- 
sorption from such high energy levels. This suggests that 
the kinetic temperature is high, in agreement with the ob- 
servations of inversion lines of ammonia up to (J, i-C )=10,10 
(which lies a t abou t 1000 K above the ground state) by 
Hcnkcl et al. (2008). On the basis of their estimates, the 
kinetic temperature Tki n is about 80 K for 80-90% of the 
ammonia gas, while the remaining fraction could even reach 
T fcm >600 K. 

A robust fitting of the CH 3 CN spectra is difficult as the 
different if-ladder transitions are close in frequency relative 
to the width of the lines. Nevertheless, a good match of 
the observations could be found using our LTE model by 
assuming a CH3CN column density of 1.3 xlO 13 cm~ 2 and 
a kinetic temperature of 50 K (see Fig lB.5p . This latter 
parameter value is co mparable to, a lthou gh slightly lower 
than that obtained bv lHenkel et al.l (|2008l ). We neglect the 
broadening of the line due to hyperfine splitting, as it is at 
most of 1.5 MHz for the principal hyperfine components of 
the Jk=§a-^a transition. 



B.14. CH 3 CN 

Methyl cyanide w as first detected in s pace toward the 
Galactic center bv ISolomon et al~ (ll97ll) . It is commonly 
detected in hot cores (e.g. lOlmi et al.ll993|) , where the heat- 
ing from massive stars is thought to evaporate the products 
of grain surface chemistry. The abundance of CH3CN in 
these regio ns of active star forma tion exceeds 10 -9 relative 
to H 2 (e.g. iKalenskii et al.ll2000D. We note that CHcC N is 
also observed in dark clouds dMatthews fc Searslll983l) an d 
has been detected in NGC253~ i|Mauersberger et al.lll99lh . 
On the other hand, CH3CN is detected in neither diffus e 
(|Liszt fc Lucasll200lD nor translucent clouds (jTurnerll200"ol) . 
with upper limits clearly lower than abundances in dark 
clouds. 

Symmetric-top molecule are excellent temperature 
probes, because radiative transitions between the differ- 
ent -fT-ladders are forbidden, so that only collisional exci- 
tation is effective. Consequently, the rotation temperature 
of the molecule can be derived within the same isT-ladder 
on the one hand, and the kinetic temperature is reflected 
in the population of the K levels on the other hand (e.g. 
IChurchwell fc HollisllT98l . 

In our LTE synthetic spectrum approach, 
the term exp(|r-^p — ) in EqJ5] was changed to 



exp 



k B T r , 



exp 



En 



k B T k 



where E 



KK 



is the low- 



est energy level of each if -ladder, to evalu ate the opacity 
of each transition (see e.g.. lQlmi et al.lfl993fl . The partition 
function was calculated accordingly by taking into account 
both the rotation and kinetic temperatures 



Q 



2_, 9I9K9J exp 

JK 



-(Ei — Ekk) 
ksTrot 



exp 



-^Y(B.l) 

KB -L kin 1 



where gj=(2J+l) and the product gigK=2 for K=0 and 
K^3n (n>0), and gigK—A for K—3n. 

Two groups of transitions of methyl cyanide fall in the 
frequency coverage of our survey, the Jk=4k-3k redshiftcd 
near 39 GHz, and the Jk=5k-^k , redshifted to ^48.8 GHz. 
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Fig. B.5. Spectra of the Ak — 3k and 5k — transitions 
of CH3CN. Our LTE prediction assuming a kinetic temper- 
ature of 50 K is overlaid. 



B.15. H 2 CCO 



Ketene was first detected in Sgr B2bv lTurnen(|l977D .It was 
detected in line surveys of regions of low-mass and mas- 
sive s tar formation (e.g . Turner! Il99lt Ivan Dishoeck et all 
119951 : iRuiterkamp et al. 20071). as well as in both translu - 



cent and dark clouds (jTurner et al.lll999l;lKaifu et al.ll2004D . 
suggesting that it is a common interstellar species. However, 
only few specific studies of the molecule have been con- 
ducted, and its formation mechanism(s) is subject to con- 
troversy: H2CCO could be either produced by gas-phase re- 
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actions (jTurner et alJ ll999) or formed b y means of surface 
chem istry, followed by grain desorption fRu iterkamp et aD 
12007ft . Its typical abundance is of the order of 10~ 10 -10~ 9 
relative to H2. 

We detect four different (ortho) transitions of H2CCO 
in our survey, yielding a rotation temperature of 2.6 ^ ' 8 K, 
significantly lower than Tcmb=5.14 K at z=0.89. 

B.16. HNCO 

Isocyanic acid was first discovered toward Sgr B2 
(ISnvder fc Buhllll972ft . We detect two lines of HNCO, the 
3o3-2o2 and the 4o4-3o3, the latter in a low- noise region of 
the survey. They do not appear to be at LTE equilibrium, 
and we could not derive a positive rotation temperature. 
The centroid velocity of HNCO is also significantly offset 
with respect to the average velocity of all molecules. We 
could not identify lines from other molecules close in fre- 
quency, and therefore discard the possibility of a false iden- 
tification. 



B.17. SiO 

Silicon monoxide is another commonly observed interstellar 
molecule, considered as a good shock tracer. Its abundance 
is low in cold dark clouds (of the order of or less than 10~ 12 
relative to H2), where it is presumably heavily depleted 
onto dust grains, but can reach values several orders of 
magnitude higher in outflows and shocked gas. It is found to 
have an intermediate a bundance (SiO/H 2 ~10~ 10 ) in diffus e 
and translucent clouds (jLucas fc Liszdl2000atlTurnerll2000h . 
Toward FG0.89SW, we find a SiO abundance of nearly 
4x 10~ 10 that is comparable to, though slightly higher than, 
that in diffuse / tr anslucent clo u ds . 

Interestingly, IZiurvs et al.l (|l989t ) found a correlation 
between ln[SiO/HCN] and the inverse of the kinetic tem- 
perature in Galactic molecular clouds: were this correlation 
to hold for the gas in FG0.89SW, the SiO/HCN ratio would 
then point to a kinetic temperature of ~30 K, a value not so 
far from the kinetic temperature estimated from CH3CN. 

The 29 SiO isotopologue is also detected and we derive 
a 28 Si/ 29 Si ratio of which is roughly half the Solar 

System value. 



B.18. CH 3 CHO 

Acetaldehyde is an asymmetric rotor, with two symmetry 
states (A and E) resulting from the rotation of the methyl 
group. Electric dipole transitions are not allowed between 
these states, which have the same weight statistic. Each 
transition is dcdoublcd by the symmetry, and the resulting 
spectrum consists of a series of closeby A/E doublets. 

The detecti on of acetaldehyd e was first reported to- 
ward Sgr B2 (|Ball et al.l Il971l) and the molecule has 
been observed in cold d ark clouds dMatthews et al.lll985l ) 
and translucent cloud s dTurner et al.l Il999f) . Interestingly, 
iChengalur fc Kanekarl ((20031) observed that CH 3 CHO is 
more widespread toward the Galactic center than most 
other complex organic molecules. They argue that the 
molecule could be produced by shock-induced disruption 
on dust mantles. 

No extragalactic detection has yet been reported for this 
species. 



B.19. NS 

Nitrogen sulfide w as first detected in Sgr B2 by 
iGottlieb etaLl (ll975h. It is detected in quiescent clouds 
(iMcGonagle et all 1994f) and in massive-star forming re- 
gions (jMcGonagle fc Ir vine 1997]). There are no reports of 
NS detection in diffuse or translucent clouds to date. 

The rotational levels of the molecule are divided into 
two levels (e/ f) with opposite parity by A-doubling, each of 
them subdivided into hyperfine structure corresponding to 
the nitrogen nucleus. We clearly detect the hyperfine com- 
ponents with highest relative intensities (see Fig |B.6[) . We 
note that the frequency interval between ^36 and 40 GHz 
is the region observed with the lowest noise level in our sur- 
vey, allowing these weak lines of NS, H 2 CS, C2S, and SO + 
to be detected. 

Taking the CS column density estimated by 
iMuller et al1(j2006l) . we find a ratio NS/ CS~0.02, compara- 
ble to the values 0.02-0.05 measured bvlHatchell. J. fc Vitil 
(|2002l ) in a sample of hot cores. iMarti'n et al.l (|2003l) 
estimated a ratio NS/CS=0.4 in NGC253, which they 
interpreted as the signature of NS chemical enhancement 
by shocks (see also lViti et al.ll2001l ). 




Sky Frequency (GHz) 

Fig. B.6. Spectrum of NS (2-1). The hyperfine structure 
with relative intensities is indicated. 



B.20. H 2 CS 

Thioform aldehyde was first de tected in absorption toward 
Sgr B2 bv lSinclair et al.l (|l973h . It is observed in star form- 
ing regions and cold dark clouds (e.g. Ilrvine et alj_[ l989: 

merl 



Min h et ai1ll991l) as well as in translucent clouds ((Turner 
119960. Tentative extragalactic detections are reported by 
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iHeikkila et all (|l999t) in the LMC and lMartm et all (|2006h 
in NGC253. 

Two lines of H2CS are detected, one of each ortho and 
para form, from which we derive an ortho/para ratio of 
3.01q 5. However, we note that another ortho line, the 2i2- 
ln at a rest frequency of 67.653820 GHz, should have been 
detected at the same intensity as the 2n-lio line, accord- 
ing to the LTE prediction. It is fitted with an opacity twice 
lower, although with a significance of just 3 a on its inte- 
grated opacity, which is the reason we discarded it from our 
estimate of the ortho-H2CS column density. 

The H 2 CO/H 2 CS ratio, derived from their para forms, 
is 32±6, fully consistent with t he present oxyg en-to-sulphur 
cosmic abundance ratio (33:1. lLoddersI 120031) . 

B.21. SO 



Sulfur monoxide was first detected by iGottlieb fc Balll 
(|1973f ) toward several s tar forming reg i ons a nd later ob- 
served in dark clouds bv lRvdbeck et al.l (|1980|) . 

From the Jr-=2i-1o line detected in our survey, we de- 
rive an abundance of ^10~ 9 relative to H2 , similar to tha t 
observed in Galactic diffuse clouds ([Lucas fc Liszt|l200l. 
and_~25 lower than that in translucent clouds (| Turner! 
Il995bh . 

The non-dctcction of the SO2 ln-Ooo ground state tran- 
sition yields a ratio SO/S02>8. Both species are observed 
in quiescent tr anslucent clou ds with an abundance ratio 
SO/S0 2 ~l-15 llTurnerlll995bh . It has been found that S0 2 
is not detecte d in diffuse clouds, do wn to a fractional abun- 
dance <10~ 9 (|Lucas fc Lisztll2002ft . 

The non-detection of 34 SO points to a n upper limit of 
SO/ 34 SO>13.7 (3ct). iMuller et al.l ([2006D derived a ratio 
32 S/ 34 S of 10.5 ± 0.6 from measurements of CS and H 2 S, 
and their isotopologues. 



B.22. SO+ 

The reactive radical SO + is characterized by a 2 n 1/ / 2 ground 
state. Its spin-doubled rotational transitions are split by A- 
type doubling. 

The SO+ J=3/2-l/2 e line is partly blended with the 
F=5/2-3/2 / hyperfine transition of NS. The width of the 
spectral feature at 1^,^36.807 GHz, as obtained by the 
fit of a single Gaussian component, is 37 ± 5 kms -1 , i.e. 
much larger than the width of other lines. Nevertheless, 
after taking into account the SO+ e/f doublet together 
with the NS line, the LTE synthetic spectrum reproduces 
fairly well the observations. 

It has been found that SO + is a rather widespr ead and 
quite abundant interstellar species ([Turner! |1994l). and i s 
seen, although weakly, in translucent clouds (|Turnerl fl996) . 
There, the abundance ratio SO/SO + is of the order of 30. 
We estimate a ratio SO/SO + =15 ± 3 of comparable order 
of magnitude toward FG0.89SW. 



B.23. C 4 H 



The b utadiynyl radical was first reported by iGuelin et al.l 
(|l978t ) toward IRC+10216 and h as been observed toward 
dark clouds bv llrvine et al.l (|1981f h It is present in translu- 
cent clouds with a quite high abundance of a few 10 -8 rel- 



ative to H2 (jTurnerl [200 01 . No extragalactic detection has 
been reported so far. 

The interaction of the unpaired electron with the molec- 
ular rotation produces fine structure, which is divided fur- 
ther into (here unresolved) hyperfine structure by the pro- 
ton spin. The resulting spectrum of C4H consists of a se- 
ries of closeby doublets, of which we detect four (7V=6-5 to 
7V=9-8) within the limits of our frequency coverage. 

B.24. HQN 

The highly unsatur ated cyanoace tylene molecule was first 
detected in Sgr B2 (|Turnerj["l971h and is observed in vari- 
ous sources, such as translucent clouds, dense dark clouds, 
and star forming regions, with abundances of ^10~ 10 -10~ 9 
relative to H2. HC3N emission has also been detecte d in 
several extragalactic objects fsee lLindberg et aT1l201ll for a 
recent review). 

We detect all four transitions of HC3N within the fre- 
quency range of our survey, from J—7-6 to J=10-9. The 
J=3-2 and J— 5-4 li nes of HC3N were pr eviously detected 
with the VLA by iMenten et all (|l999h . More recently, 
iHenkel et al.l (|2009f ) reported the detection of all transitions 
from J=3-2 to J— 7-6 with the Effelsberg 100 m telescope. 
The J=10-9 line was not detected. The results of their ob- 
servations are reported in the rotation diagram together 
with our data points (Fig[3]). They appear to be incon- 
sistent with one another. The integrated opacities of the 
J=7-6 line, observed in common with both works, differ 
by a factor of ~2, as for their fitted linewidths. Time vari- 
ablity and/or errors in the spectral baseline could explain 
the observed differences. We note that the higher S/N J=7- 
6 line appears slightly asymmetrical, with a shoulder on 
the blue-shifted side, reminescent of that of t he HCO + 2-1 
profile observed at higher sp ectral resolution (|Muller et al.l 
120061; IMuller fc Gu clin 200j|. This transition was therefore 
not included in the Gaussian fit to derive the velocity and 
linewidth. 

The points corresponding to the four lines observed by 
us appear remarkably well-aligned in the rotation diagram, 
suggesting that the lines are at LTE. They yield a rotation 
temperature of T rot =6.3 ± 1.3 K, just slightly higher than 
the expected CMB temperature. On the ot her hand, the 
points corresponding to the lines observed bv lHenkel et al.l 
(|2009f ) at lower frequencies are not aligned in the rotation 
diagram. While we cannot discard possible time variations, 
another explanation could be a change in th e source cov- 
ering factor with frequency, as discussed by IHenkel et alj 
(|2009f ). Observations of all these transitions within a short 
period of time would help clarify the interpretation of these 
data. 



B.25. C 2 S 

The thioxo e theny lidene radical was first identified by 
ISaito et all ([l987) in previous spectra of Sgr B2 and 
TMC 1, based on laboratory spectroscopy studies. It has 
been found that C2S is present in translucent clouds 
([Turner et al.lll998h with a typical abundance of ^10~ 9 and 
i s dete cted in the 2 mm survey of NGC253 bv lMartm et~aT1 
(2006). It is the heaviest species detected in our survey. We 
measure an abundance of ^2xl0 -10 relative to H2. 
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Few lines (10 > J u > 16) of C 3 S fall within the fre- 
quency coverage of our survey. Their non-detection gives 
a (not tightly constraining) ratio C2S/C3S>0.7. Search 
for lower J t ransit ions of C3S would be more sensitive. 
IFuente et al.l (|1990f ) measured a typical abundance ratio 
C2S/C3S=5 in dark clouds. 



